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ABSTRACT 


RAPLOT II is a computer program for processing radiation and naviga- 
tion data from field surveys of the Radioisotopic Sand Tracer (RIST) study, 
but is applicable to any survey type operation on the nearshore Continental 
Shelf. Collected data are punched onto paper tape by the data collection 
computer on the research vessel. The data are later transferred to mag- 
netic tape which provides the input for the RAPLOT II Program. Program 
control parameters are on punched cards. The navigation data, which con- 
sists of ranges to two shore-based radar beacons, are first edited for 
spurious data, and then converted to rectangular coordinates (in this case 
the California Lambert Coordinate System). Radiation data are converted 
to count rate as counts per second. Background count rate is computed 
and subtracted from the observed count rate, and any radiation counts that 
are significantly above the background count rate are corrected for time 
of decay since the isotope was injected. Output from the program is in 
three forms - printed output, graphical output, and magnetic tape record. 
The processed data are transferred to magnetic tape and made available for 
further processing such as the generation of contour maps. 


FOREWORD 


CERC Miscellaneous Paper 2-69, Radioisotopic Sand Tracer Study, Point 
Conception, California, reported the early results of the RIST study. This 
study is part of Contract AT(49-11)-2988 between the Atomic Energy Com- 
mission and CERC. Other participants in this continuing multi-agency study 
are the Oak Ridge National Laboratories of the Atomic Energy Commission; 

U. S. Navy Pacific Missile Range; U. S. Air Force (Western Test Range, First 
Strategic Aerospace Division); U. S. Army Corps of Engineers Los Angeles 
District office; NASA (Nuclear Systems and Space Power Division), the State 
of California (Department of Water Resources) and U. S. Army Mobility Equip- 
ment Command. The study involves the collection and analysis of an enormous 
amount of data. Computer processing is the only means by which these data 
could be handled. 


Philip A. Turner, a geologist, developed the original RAPLOT program 
and prepared this report. The work was done under the general supervision 
of David B. Duane, Chief, Geology Branch, and George M. Watts, Chief, En- 
gineering Development Division. CERC continues to refine data processing 
and improve the printed and graphic output which consists of maps drawn by 
an incremental plotter. 


At the time of publication, Lieutenant Colonel Edward M. Willis was 
the Director of CERC; Joseph M. Caldwell was Technical Director. 


NOTE: Comments on this publication are invited. Discussion will be 
published in the next issue of the CERC Bulletin. 


This report is published under authority of Public Law 166, 79th 
Congress, approved July 31, 1945, as supplemented by Public Law 172, 88th 
Congress, approved November 7, 1963. 
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Section A. INTRODUCTION 


In 1966 the Coastal Engineering Research Center (CERC), in coopera- 
tion with the Atomic Energy Commission, initiated a 3-year radioisotopic 
sand tracer study of littoral transport around Point Conception, Cali- 
fornia. The purpose was to develop and use radioactive tracers for 
research in sand movement and littoral processes. The objectives 
included determination of suitable radioactive isotopes, development of 
handling and survey procedures, and development of computer programs for 
editing, processing and graphical display of the data. At the same time, 
studies of sediment transport around the Point Conception headland and 
of the mechanics of littoral transport were conducted. Methods developed 
by this program have direct application to engineering design of harbor 
development and beach erosion prevention, and quasi-military application 
such as the location of radioactive and other toxic materials. 


Sand grains indigenous to the study area are labeled with a radio- 
isotope that does not adversely affect their hydraulic properties. A 
mobile detector system, using cesium iodide crystals and housed in a 
"ball" towed behind an amphibious vehicle, detects the location and 
intensity of the radiation. Concurrently, additional field data are 
collected on sediment size and composition, isotope distribution, beach 
and nearshore bottom topography, weather, and sea and swell conditions. 


During a sand-tracer field investigation, radiation measurements are 
made continuously as the mobile detector system is towed along a beach, 
through the surf, and over the offshore bottom. With a time selection 
mode for data acquisition available in increments from 0.1 to 10.0 
seconds, a large mass of data is accumulated in a few hours. During a 
field test, surveying may go on several hours a day for several weeks. 
Computer processing is necessary to study and evaluate the great volume 
of collected data. Plotting and posting of the survey data is also use- 
ful for monitoring field operation. 


The initial field investigations at Surf, California, relied on 
manual preparation of maps and subjective interpretation of data printed 
by the teletype of the onboard data acquisition system. It immediately 
became evident that computer processing and plotting must be employed in 
future operations, and ‘CERC undertook development of computer programs 
to generate plots useful to continuing field operations. This specific 
program is called RAPLOT. The first version of the program was used to 
process the data collected in the December 1967 field test at Point 
Conception, California, and is documented in Appendix D of CERC report, 
M.P. 2-69, by Duane and Judge (1969). RAPLOT was originally compiled 
and made operational on the UNIVAC 1108 at the National Bureau of Stand- 
ards, with graphic display on a Benson-Lehner incremental plotter at CERC. 
Later, the program was adapted to the IBM 7094 Computer and Milgo plotter 
at the Western Test Range of Vandenberg Air Force Base. This version of 
the program was used to support later field tests at Point Conception 
and Surf, California, in September and October 1968 and in February and 
June 1969. 


Improvement and modification of program data acquisition systems 
created changes in the collection format which required a new version 
of the RAPLOT program. At the same time, experience gained from running 
the RAPLOT program was incorporated to provide an improved and more 
sophisticated data processing capability for the RIST project. The new 
program is called RAPLOT II and was used at CERC for processing field 
data for all tests after December 1968. Subsequently, the program was 
modified to run on an IBM 7094 computer to produce plots on a Stromberg- 
Carlson 4020 cathode ray tube. This version is called RAPLOT III (See 
Appendix B, page 47). 


Section B. SUMMARY OF RAPLOT II PROGRAM 
1. Hardware Requirements 


RAPLOT II was written in FORTRAN V for the UNIVAC 1108 Computer 
and EXEC II operating system at the National Bureau of Standards (NBS) 
in Gaithersburg, Maryland. CERC is connected to NBS by a leased tele- 
phone line and has a UNIVAC 1004 for its remote terminal. The NBS 
UNIVAC 1108 has 65,000 words of core memory of which 38,400 are re- 
quired for the RAPLOT II Program. FORTRAN V employs advanced features 
not found in standard FORTRAN IV. They are the NTRAN subroutine for 
executing binary input-output commands, and the FLD function, a bit 
manipulation routine. 


For field program use, RAPLOT II has been modified and written in 
FORTRAN IV to run on an IBM 7094 computer at the Pacific Missile Range 
data processing center at Point Mugu, California. The size of the 
program had to be reduced considerably to fit into the 32,000 word 
memory. Consequently, the processing of the data is less thorough, and 
the writing out of the processed data on magnetic tape was eliminated. 
This version was used to support a RIST field test near Point Mugu where 
the primary requirement was for quick printout and graphic display of 
the processed survey data. 


2. Program Description and Logic 


A flow chart of the program is given in Figure 1; an outline 
description of the main steps in the program follows: 


a. Read in two program control cards and a file legend card. 
Additional control parameters are computed from these input parameters. 


b. Read in from magnetic tape a data file from a radio- 
isotopic tracer survey. 


c. Check radar beacon ranges for errors. If distance-time 
ratio for successive ranges indicates a ship speed greater than 6 knots, 
or 3.09 meters per second, correct the ranges by linear interpolation. 
Experience has indicated that these beacon ranges may be erroneous 
(as much as an order of magnitude) as often as 5 percent of the time. 
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d. Compute the position of the survey vessel at each fix. The 
program control cards provide the California Lambert Coordinates of the 
shore beacon, and the data input tapes provide the beacon ranges. By the 
cosine law this can be translated into distance in terms of a rectangular 
coordinate system: 


= (DB2 + D* - D2)/(2°DB) 
DX, (DB* + DF 5) 
rane 2) 35 
= (Dab X 2 
DY, ( i D 
Where Dj, = the distance of survey vessel from the upcoast radar beacon 
D, = the distance from the downcoast beacon 


DB = the distance between the two beacons 


DX, = the distance of the fix position of survey vessel from 
the upcoast beacon in a direction parallel to a line 
intersecting the two beacons 


DY] = the distance of the fix from an upcoast beacon in a 
direction normal to a line intersecting the two beacons. 


These distances may be translated and rotated to give the Lambert 
Coordinates for each fix 


NORTH ¢ = DX) sin6é + DY, cos§ + NORTH, 


EAST ¢ = DX} cos6é - DY} sind + EAST, 


Where 6 = the angle of rotation of the coordinate system 
NORTH, and EAST, = the Lambert coordinates of the upcoast 
radar beacon 
NORTH, and EAST, = the coordinates of the fix. 


Occasionally the two radar beacon ranges will produce an imaginary 
position. When this occurs, one of the coordinates of the radar fix is 
defined by the square root of a negative number. When such an imaginary 
fix occurs, the line of data producing the imaginary fix is eliminated 
from the file being processed. 


e. The north and east coordinate arrays are searched for ex- 
treme values by Chauvenet's criterion and any data set containing extreme 
values is eliminated. Chauvenet's criterion is based on the normal dis- 
tribution, and a value is rejected if the probability of occurrence of 
such a deviation from the mean of the n measurements is less than 1/2 n 
(Po Wo Relekere5 IK) © 


f. Plot the track of the survey vessel by a call to subroutine 
BENLH2. 


g. Compute the position of the mobile underwater detector for 
each fix by correcting for the distance the detector is towed astern of 
the survey vessel. Cable length, mean water depth, and length of survey 
vessel, are all input on the program control card. Assuming a position 
for the detector vehicle at the beginning of the survey, the position of 
the detector for each fix is computed by linear interpolation from the 
present position of the vessel to the last interpolated position of the 
detector. 
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where XD and YD are the coordinates of the detector and XV and YV are 
coordinates of the vessel. CABLE is the horizontal distance from the 
detector to the survey vessel. 


h. Radiation data are converted to counts per second, and the 
mean and standard deviation of the background count rate is computed. 
An estimated background count rate is entered on the data control card 
for each channel. This estimated background count rate is used to 
compute the extreme values of the range of the background radiation 
level, again by means of Chauvenet's criterion. All counts between 
these extreme ranges are averaged to obtain the mean background count 
rate for each channel. 


i. Plot uncorrected radiation values by a call to BENLH2 
through RDPLT2 entry. 


j. Correct the radiation data by subtracting the mean back- 
ground count rate. Then correct the remainder, if it is significantly 
higher than background, for the time of the decay since the isotope was 
injected. 


CCR = (NCR - 3c)e4t + 30 


loge (2) 
where he ie D and 
2 
where CCR = the corrected count rate 


NCR = the net (observed radiation value less mean background) 
count rate 


o = the standard deviation of the background count rate 
t = the elapsed time from the injection to the time of 
the fix 
T, = the isotope halflife in hours. 


1 
2 


Only radiation count rates that are significantly greater than back- 
ground are corrected for time of decay. Otherwise, in a survey made one 
or more halflives after the injection, the correction would inflate back- 
ground radiation readings to such a degree that they would appear to be 
significant. 


k. Compute the weighted mean and weighted standard deviation 
of the coordinate location of the radioactivity. Compute the 95 percent 
confidence limits of the mean radiation location. 
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1. Print out the sequence number, time, radar beacon ranges, 
coordinates, and activity for each line of data. 


m. Return to the beginning of the program to read in more data 
control cards, and process another file of data unless it is the end of 
job, in which case processing ends. 


Section C. SUBROUTINES USED IN THE PROGRAM 
1. BENLH2 - Plotting Subroutine 


Subroutine BENLH2 performs the operations necessary to produce 
a plot of the trackline and of the radiation data on the Benson-Lehner 
incremental plotter. BENLH2 does this by calling the several subroutines 
of the Benson-Lehner plot package which translates the data given to the 
subroutines by BENLH2 into plot commands that are written out onto mag- 
netic tape. The tape is used to drive a Benson-Lehner model 305 incre- 
mental plotter off-line. Figure 2 is a flow chart of this subroutine. 


Subroutine TRACK is substituted for subroutine BENLH2 to produce 
the RAPLOT III program. TRACK interfaces the RAPLOT program with the 
plot subroutines for the Stromberg-Carlson 4060 cathode ray tube. Other- 
wise, it functions essentially like subroutine BENLH2. Both subroutines 
plot the trackline followed by the survey vessel, and both will also plot 
the uncorrected or corrected radiation values for each survey. These are 
symbol plots in which the count rate is represented by a symbol indicat- 
ing a value between arbitrarily fixed limits. For uncorrected radiation, 
the limits are established in terms of standard deviation from the mean 
background count rate. For corrected radiation count rate, the internal 
limits for each symbol are established on a power of 2 scale. Also, a 
special symbol is used to indicate background if the count rate is less 
than 3 standard deviations from the background count rate. If the count 
rate is more than 3 standard deviations below background, the value is 
not plotted at all. The reason for this is that an abnormally low count 
rate may indicate that the detector was ''flying"’ meaning that it was not 
in contact with the ocean bottom at the time of the fix. Both BENLH2 
and TRACK will plot the location of the mean radiation position referred 
to as RADBAR. They may also plot the location of the injection site, if 
this is desired. 


2. Statistical Subroutines 


The subroutines described below are used to perform certain sta- 
tistical computations required by the RAPLOT II program. To reduce the 
time necessary to prepare the program, these subroutines were taken from 
the STAT-PACK, a library of statistical subroutines available on the 
UNIVAC 1108 system. Further details may be obtained from the STAT-PACK 
Programmers' Reference Manual. 


Subroutine STDEV computes the standard deviation of an array by 
the following formula: 
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Figure 2. Flow Chart of BENLH2 Program 


where X = the array of values 


>| 
u 


the arithmetic mean of the X array 
n = the number of elements in the X array 
o = the standard deviation 


Subroutine AMEAN is called to compute the arithmetic mean of the X 
array. 


Subroutine TINORM computes the value of the inverse of the normal 
distribution by a rational approximation. The inverse normal distribution 
is defined by the solution for X of the following integral equation 


xX 2, 
e are dt =a 


1 
a 


where a is the probability for which the ordinate is to be calculated. 
The rational approximation itself may be ascertained from the FORTRAN 
statements in the subroutine. 


3. Input-Output Subroutines 


The binary input-output statements in FORTRAN, although conven- 
ient to use, cannot make efficient use of magnetic tapes because the 
language does not permit total parallel processing. Furthermore, a 
considerable amount of time is used in processing a list because of its 
generality. More efficient use is attained by buffering the input or 
output in parallel with computing. On the UNIVAC 1108 system the NTRAN 
subroutine provides a means of buffering through a call statement in 
FORTRAN: 


CALL NTRAN (UNIT, sequence of operations) 


in which UNIT is an integer constant or variable designating the logical 
unit. If the unit is not busy, NTRAN initiates the first operation and 
stacks the remainder in a waiting list. If the unit is already busy, 
then the entire sequence is stacked in a waiting list and chained to the 
previously stacked sequence. Control is then returned to the program 
following the call NTRAN statement. When an interrupt occurs, NTRAN 
records the transmission status, initiates the next operation in the 
chain and returns control to the interrupted program. 


The NTRAN subroutine is used to write out on magnetic tape the 
arrays containing the legend, sequence number, underwater detector co- 
ordinates and corrected count rate for both radiation channels for each 
file of survey data being processed. The information thus stored on 


magnetic tape may be used for further processing such as drawing contour 
maps. No provision is made for outputting the processed data in RAPLOT 

III because no further processing was planned for any of the field test 

sites. Also deletion of the output statements reduces the running time 

of the program. 


Call NTRAN statements are also used in RAPLOT II to position the 
input tape by end of file marks. 


Section D. PROGRAM INPUT 


I, Geaeal Input 


Input for the RAPLOT II program comes from punched cards and 
from magnetic tape. Three program control cards are read in for each 
file of survey data processed. Formats, and descriptions of the vari- 
ables written onto the cards are given in Table I; names given to the 
variables in this table are the ones employed by the program. Figure 3 
shows a data sheet that is used for filling out the control cards. It is 
useful not only for filling in the parameters on the cards, but also for 
keeping track of the data files when large numbers of files are being 
processed. 


The last 12 characters of the legend (Card 3, spaces 66-78) are used 
for file identification. Before a file of survey data is read in, a 7- 
word identification array called SENTNL is read in from the beginning of 
the tape file. The first two words in SENTNL are compared with the last 
two words in LEGEND. If a match is not found, the tape is positioned at 
the beginning of the next file and a new SENTNL array is read in. This 
procedure is in lines 48 through 52 of the source language listing of 
RAPLOT II in Appendix A. It has been deleted from RAPLOT III. 


Be Tape Input 


When a RIST survey is underway in the field, data from various 
sensors are assembled by the onboard detector system, and punched on 8- 
channel paper tape in American Standard Code for Information Interchange 
(ASCII) code. At presént, there are seven data fields for each line of 
record. A brief description of the variables and the tape record format 
is given in Table II; Figure 4 is a sample listing from a paper-tape 
data file. 


In the actual processing of thé RIST data files, it has been found 
necessary to edit the data files prior to putting them through the RAPLOT 
II program. The edited data are written out in unformatted magnetic tape 
files. This is why the READ statement in line 58 of the listing in 
Appendix A is an unformatted FORTRAN READ statement. Unformatted (or 
binary) input-output statements are much more efficient for the computer 
to execute than formatted statements. For that reason, they are employed 
wherever possible. RAPLOT III employs a formatted READ statement (line 
51, Appendix B) for inputting a data file, since the need for short turn- 
around time is greater than the need to edit data files when supporting a 
field test. 


TABLE I 
Format of Program Control Cards 


I Data Control Card (3F3.0,2F7.0,1X,2F2.0,F3.0,F5.2,3F2.0,4(1X,F7.0),13) 


Column Variable Description 
1-3 CABLE Length of cable, in feet, on which the detector 
is towed. 
4-6 DEPTH Average water depth plus freeboard, in feet. 
7-9 BOAT Distance from radar mast to cable stanchion 
in feet. 
10-16 BKG(1) Estimated background count rate (counts per 


second) for radiation channel 1. 
17-23 BKG (2) Same for radiation channel 2. 


25-28 ZHR, ZMIN The time of injection in hours and minutes 
(24-hour clock). 


29-31 DAYS The number of days that have elapsed since 
the injection. 
32-36 HLIFE The halflife of the radioisotope in days. 
37-40 SETIME, 
RMIN The time when the survey was started, in hours 


and minutes (24-hour clock). 


41-42 SEC The time, in seconds, between successive fixes. 
44-50 BEACIN 

52-58 BEACIE California Lambert Coordinates of the radar 
60-66 BEAC2N beacons. BEACON 1 is always the upcoast beacon. 
68-74 BEAC2E 

75-77 ISKIP The number of lines to skip at the beginning of 


a data set in order to avoid reading in some 
x bad data. 


TABLE I (Continued) 


Format of Program Control Cards 


II Plot Control Card (411,F10.0,1X,312,1X,F10.0,2(1X,F10.0),1X,A6) 


Column Variable Description 
1-4 NOPT Plot option controls. A numeral 1 in the 


column indicated causes the various options 
to be executed. 


1 Plot trackline followed by survey vessel. 

2 Plot uncorrected radiation values. 

3 Plot radiation values corrected for background 
and for decay since injection. 

4 Unused, 

5-14 SCALE Map scale in feet per inch. 

16-17 NPLT (1) Options used for spotting data points for each 

18-19 NPLT (2) plot option. User can specify that every Nth 

20-21 NPLT(3) point be plotted. If left blank, every point 


will be plotted. 


23-32 GRID Intervals on the coordinate grid at which 
Lambert Coordinates will be posted. If 
GRID=0, tick marks are not plotted. 


34-43 DUMP 1 Distance in meters from the injection site to 

45-54 DUMP 2 the upcoast and downcoast beacon, respectively. 
If the fields are blank, then the injection 
site is not plotted. 

56-61 INDATE Day, month and year that the radioactive sand 


was injected. 
III Plot legend card (13A6,L2) 


1-78 LEGEND A descriptive legend that is included on the 
printed output, and is also written on the 
lower margin of the map. 


79-80 JOBEND The letter 'T' is entered here if the data file 
being processed is the last one in the job. 
Otherwise the field is left blank. 
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CARD 2 PLOT CONTROLS 


CARD 3 LEGEND 


Data Sheet for RAPLOT II Control Cards 


Figure 3. 


TABLE II 


Format of Paper-tape Record on Which RIST Survey Data is Collected 


Paper-Tape Format (16,3F7.1,3F7.0) 


Column Array Description 

1-6 NMBR A line (record) number which is incremented 
with each line of data unless manually reset. 

8-13 TIME Cumulative time, in tenths of seconds, since 
the beginning of the survey. Time is reset 
to zero at the beginning of each survey. 

15-20 Dy Distance, in tenths of meters, from the 
upcoast radar responder beacon. 

22-27 D, Distance, in tenths of meters, from the 
downcoast radar responder beacon. 

29-34 RAD Radiation Channel 1: pulses from differential 

1 : ae : 

discriminator accumulated for time between 
successive records. It is automatically 
reset to zero at the writing of each record. 

36-41 RAD, Radiation Channel 2: same description. 

43-48 FATH Depth of water recorded by fathometer. This 


is not presently in operation. 


000061 001240 001763 004900 000236 000201 000000 
000062 001260 001793 004916 000175 000142 000000 
000063 001280 001814 004919 000155 000132 000000 
000064 001300 001846 004924 000150 000130 000000 
000065 001320 001861 004929 000142 000119 000000 
000066 001340 001878 004930 000152 000136 000000 
000067 001360 001920 004920 000149 000124 000000 
000068 001380 001925 004825 000152 000116 000000 
000069 001400 001958 004947 000136 000114 000000 
000070 001420 001980 004975 000130 000110 000000 
000071 001440 001022 004932 000130 000101 000000 
000072 001460 002027 004932 000127 000104 000000 
000073 001480 002043 004945 000144 000117 000000 
000074 001500 002092 004977 000145 000123 000000 
000075 001520 002120 004966 000108 000089 000000 
000076 001540 002143 004970 000140 000115 000000 
000077 001560 002157 004970 000143 000117 000000 
000078 001580 002189 004979 000174 000136 000000 
000079 001600 002206 004981 000142 000114 000000 
000080 001620 002234 004999 000138 000122 000000 
000081 001640 002265 004996 000133 000108 000000 
000082 001660 002286 004998 000188 000155 000000 
000083 001680 002315 004016 000122 000101 000000 
000084 001700 002327 005048 000164 000131 000000 
000085 001720 002365 005011 000166 000133 000000 
000086 001740 002386 005021 000178 000151 000000 
000087 001760 002403 005039 000152 000120 000000 
000088 001780 002448 005075 000133 000105 000000 
000089 001800 002471 005073 000154 000131 000000 
000090 001820 002497 005077 000169 000143 000000 
000091 001840 002509 005089 000151 000128 000000 
000092 001860 002549 005109 000174 000153 900900 
000093 001880 002580 005135 0090153 000130 000000 
000094 001900 002611 005138 000158 000126 000000 
000095 001920 002629 005156 000154 000124 000000 
000096 001940 002658 005164 000145 000114 000000 
000097 001960 002693 005191 000170 000141 000000 
000098 001980 002719 005207 000189 000151 000000 
000099 002000 002757 005211 000167 000137 000000 
000100 002020 000774 005219 000147 000120 000000 
000101 002040 000808 005238 000170 000133 000000 
000102 002060 002836 005266 000173 000129 000000 
000103 002080 002868 005272 000149 000126 000000 
000104 002100 002902 005280 000146 000126 000000 
000105 002120 002916 005302 000192 999149 000000 
000106 002140 002958 005322 000140 000114 000000 
000107 002160 002997 005342 000170 000136 000000 
000108 002180 002028 005355 000163 000133 000000 
000109 002200 003044 005363 000142 000121 000000 


Figure 4. Sample Listing of Input Data File for RAPLOT I1 
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Section E. PROGRAM OUTPUT 


1. Printed Output 


Three forms of output are generated by the RAPLOT II program - 
printed output, graphical output, and magnetic tape output. 


The printed output for one data file consists of one page listing 
the program control parameters and summary statistics (Figure 5) followed 
by a complete listing of the unprocessed and processed data arrays 
(Figure 6). 


The printout of program parameters and summary statistics for RAPLOT 
III (Figure 7) is virtually the same. However, the radar-range data and 
uncorrected radiation data are deleted from the printed listing (Figure 8). 
Also, only one channel of corrected radiation data is included. 


Because of the large volume of data files that may be processed by 
the RAPLOT II program, it is usually necessary to keep track of the number 
of pages of output. Fifty lines are printed to a page of output; one page 
is required for program control parameters and another page for summary 
Statistics. Therefore, the number of pages per data file = 2 + (nunber of 
lines of data/50.) 


2. Graphical Output 


Graphical output from RAPLOT II is in the form of plot commands 
that are either written out on magnetic tape or punched onto cards to be 
used to drive a Benson-Lehner plotter off-line. Under the present setup, 
the plot commands are written out on logical tape unit 9; therefore, 
this same unit cannot be used for other tape input or output without first 
making changes. On the UNIVAC 1108 system, the plot commands are blocked 
and written out in card image length records at 556 bits per inch in even 
parity, BCD. Experience has shown that one full tape will hold approxi- 
mately 20 separate plots. Assuming that 2 plots are generated per file 
of data, this means that 10 files of radiation survey data will generate 
a full tape of plot commands. 


RAPLOT II graphical output consists of three types of plots - plots 
of the trackline followed by the survey vehicle, plots of uncorrected 
radiation data(for plotting background radiation), and plots of cor- 
rected radiation data. Selection of the various plots is controlled by 
the variable NOPT in Table I. The plot selection options are independent 
of each other. Any one of the three plot types may be selected, or all 
three if it is desired. Normally, the trackline plot is selected and 
then, depending on whether the data is from a background or a radiation 
survey, either the plot of uncorrected radiation or corrected radiation 
data is selected. Plots of the survey vehicle trackline and the cor- 
rected radiation data are shown in Figures 9 and 10. 
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PT MUGU GROIN SITE BACKGROUNT SURVEY 1 


LINE TIME DISTANCE TO BEACON BOAT COGROTNATES 


SEC 1 2 NORTH 

t) 2. 228177. 
1 4 228192. 
2 be 226182. 
3 8. 228181. 
4 10- 228171. 
Seize 2281686 
6 14. 228160. 
T 166 228145. 
8 18 226145. 
S$ 206 228141. 
19226 228130. 
1k 24. 2281186 
12 266 228119. 
130 (28. 228124. 
14 30. 228198. 
15e 32s 228105. 
16 34. 228191. 
18 38. 228072. 
19° 40. 228063. 
20 42.6 228050. 
21 44. 228040. 
22 466 228023. 
23 48. 228015. 
24 56 228011. 
25 52. 228995. 
26 54 227994. 
27 56.6 227989. 
28 «58. 2279926 
29 «of. 227970. 
39 62.6 2279706 
31 64 227960. 
32 666 227957. 
33 68. 227942. 
34 70. 227935. 
35 72. 2279196 
36 14. 227921. 
37 16. 227913. 
38 9078 2279146 
39° BN 227906. 
49 82 227895. 
41 Bae 2278946 
42 86e 227881. 
43 88. 227886. 
44 9M. 227875. 
45 920 227862. 
46 9b. 227854. 
47 %6o 2278560 
48 98. 227845. 
49 100. 227833. 
50 102. 227825. 


EAST 
1648434. 
1648435. 
1648427. 
1648424. 
1648419. 
1648416. 
1648408. 
1648406. 
1648492. 
1648398. 
1648394. 
16483926 
1648385. 
1648379. 
16483736 
1648361. 
1648354. 
16483356 
16483256 
1648316. 
1048308. 
1648391. 
1648295. 
1648284. 
1648275. 
1648273. 
1648264. 
16482546 
1648251. 
1648241. 
1648232. 
1646220. 
16482146 
1648294. 
1648295. 
1648199. 
1648197. 
1648192. 
1648135. 
1648179. 
16461736 
16481736 
1648164. 
1648157. 
1648152. 
1648147. 
16481416 
1648138. 
1646134. 
1648128. 
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BALL COORDINATES 

NORTH EAST 
228092. 1648429 
228196. 16484396 
228106. 1648430. 
2281065 1648430. 
227106. 1648439. 
228196. 16484306 
228196. 16484306 
228196. 16484306 
2281065. 1648430. 
22A196¢ 1648430. 
228106. 1648430. 
228196. 1648430. 
22A1N6. 16484396 
226106. 16484396 
2281065 1648430. 
228106. 16494306 
228106. 164A43N. 
22AlClL. 1648416. 
223096. 16484046 
228990. 1648392. 
228084. 16483826 
228075. 16483706 
229969. 1648362. 
229963- 16483556 
228957. 1648347. 
228950. 16483396 
228544. 16483326 
224040. 16483256 
228029¢ 1648314. 
228024. “1648307. 
228915~¢ 16462986 
228099. 1648288. 
227999. 16482766 
227992. 16482696 
227983. 1648262. 
227982. 1644260. 
227976. 1648255. 
2279742 16482536 
227967. 16482466 
227958. 1648237. 
2279550 16482340 
227947. 16482286 
2279452 16482266 
2279360 1648217. 
227927. 1648208- 
2279202 1648293. 
227918. 1648201. 
2279100 1648194. 
227901le 1648187. 
227893. 1648181. 
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RADIATION COUNTS/SEC 
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FEET 


Figure 8. Sample of Printed Output of RAPLOT III Processed Data 
for One RIST Data File 
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SCALE IN FEET 


—————— INJECTION DATE 
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Figure 9. Trackline Plot Produced on Benson-Lehner Plotter by 
Subroutine BENLH2 of RAPLOT II 
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if “Pr zp ats + 
CORRECTED COUNT RATE (C/S) 
y. OOR LT. OR EQ. “61.1 
©, COR.LT. OR .€Q. +158.7 

CCR .LT. OR .€Q. -317.5 

COR .LT. OR EQ. -635.0 

CCR .LT. OR .€Q. +1270.0 

CCR .LT. OR EQ. +2539.9 

of * o COR.LT. OR EQ. +5079.9 

FE : COR GT. 5079.9 


ema Le ae! ee 4 


+444900 
+1216800. 


fo 
16800 +1217000. +1217500. +1218000. +1218500. , 
SURF 600 FT S R-158 0/S RADIATION SURVEY AU-198 1451 10/04/1968 5 


SCALE IN FEET PLOT OF CORRECTED RADIATION/ INJECTION DATE 
eager Oe eR T3000 7400 BACKGROUND RADIATION RATIOS 041068 


Figure 10. Plot of Corrected Radiation Data Produced on Benson-Lehner 
Plotter by Subroutine BENLH2 of RAPLOT II 
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The graphical output from RAPLOT III is produced on the Stromberg- 
Carlson 4020 computer recorder. The technique for getting the plots is 
somewhat the same as the Benson-Lehner plotter in that the plot commands 
are written out on magnetic tape and are then used off-line to generate 
the plots. There is a highly developed software package that goes with 
the S-C 4020 and if it is planned to use this method of generating plots, 
the Programmers' Reference Manual for the S-C 4020 should be used. Figure 
ll is a trackline plot produced on the computer recorder at the Pacific 
Missile Range data processing center at Point Mugu. Figure 12 is a plot 
of uncorrected radiation data that was also produced there. 


3. Magnetic Tape Output 


The processed data is stored on magnetic tape for future refer- 
ence by the calls to the NTRAN subroutine. The arrays written out are: 
LEGEND, sequence number (NMBR), coordinates of radiation location (NCORD 
and ECORD), and corrected radiation data (CCR). This procedure has been 
eliminated from RAPLOT III. 


Section F. INSTRUCTIONS FOR RUNNING PROGRAM 


An example of a job deck setup for running RAPLOT II is shown in 
Figure 13. Further instructions on running jobs are in the 1108 EXEC II 
Programmers' Reference Manual. Running time depends on the number of files 
being processed and the number of records in each file. Figure 14 shows 
a graph of UNIVAC 1108 central processing unit time to process one file 
versus the number of records in the file. 


For running version III on the IBM 7094, consult the IBM reference 
manuals for FORTRAN IV and the system monitor (IBSYS). It may also be 
useful to have a copy of the Programmers' Reference Manual for the 
S-C 4020 Computer Recorder. 


Section G. RADIATION CONTOURING PROGRAM 


Present programming effort is directed toward completing RADCON, a 
FORTRAN V program for drawing contour maps of radiation data. The input 
for this program will be the files of processed data on magnetic tape that 
have been generated by RAPLOT II. If necessary, two or more data files may 
be combined to provide the input for one contour map. The radiation data 
are first smoothed by a moving average (LINAVE) procedure and then inter- 
polated over a uniform grid by a weighted least-squares numerical approxi- 
mation (NUPRX). The resulting grid is then contoured at equally spaced 
intervals of the gridded values. There are also options for transforming 
the radiation data, although the number of options and types of trans- 
formation is still undecided. 
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Figure 11, Trackline Plot Produced on S-C 4060 Computer Recorder by 
Subroutine TRACK of RAPLOT III 
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Figure 12. Plot of Uncorrected Radiation Data Produced on S-C 4060 
Computer Records by Subroutine TRACK of RAPLOT III 
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BF IN 
BFIN 
BEOF 


V 
PROGRAM CONTROL CARDS 
BN XQT RAPLOT 


WV Norw FORTRAN DECK 
BN FOR TINORM, TINORM 


Vs 
J KMEAN FORTRAN DECK 
BN FOR AMEAN, AMEAN 


Ye 
Y: 
YSTDEV FORTRAN DECK 
BN FOR STDEV, STDEV 


YL 
Ye 
V2 
YMBRENLH2 FORTRAN DECK 
BN FOR BENLH2, BENLH2 


Y 
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YRAPLOT I FORTRAN DECK 
BN FOR RAPLT2, RAPLT2 
INF PLOT (LOAD BENSON-LEHNER plotting routines ) 
8B XQT CUR 
BNKE ASG C=CERCI19 (output tape for plot commands ) 
BN ASG B=CERCI21 (output tape for stored dota ) 
BAEK ASG A=CERCI75 (INPUT TAPE ) 


$C RUN TURNER, 21053, 10, 200 


Figure 13. Example of a RAPLOT II Job Deck Setup for UNIVAC 1108 
Running under EXEC II 
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APPENDIX A 
LISTING AND INDEX OF RAPLOT II PROGRAM 
The source deck has been analyzed routine-by-routine and a directory 
of all statement numbers and variable names used in the source deck. The 
symbols are listed numerically and alphabetically. Beside each symbol 
appears line numbers followed by a code which designates where and how 
the symbol was referenced. The codes used are as follows: 


Code Meaning 


Blank A simple reference, such as the use of a statement number in an 
IF statement or a variable used in an arithmetic statement. 


= Arithmetic definition of a variable which appears on the left 
of an = sign. 


* . Statement number defined at this line. 


AG Variable appears as an argument in a subroutine or function 
statement or as an argument in a CALL statement. 


co Variable appears in a COMMON statement. 

cx Variable appears in a COMPLEX statement. 

DA Variable appears in a DATA statement. 

DB Variable appears in a DOUBLE PRECISION statement. 

DI Variable appears in a DIMENSION statement. 

EQ Variable appears in an EQUIVALENCE statement. 

EX Function name appears in an EXTERNAL statement. 

IN Variable appears in an INTEGER statement. 

LG Variable appears in a LOGICAL statement. 

NM Variable appears in a NAMELIST statement. 

PR Variable or FORMAT number appears in a PRINT statement. 
PU Variable or FORMAT number appears in a PUNCH statement. 
RD Variable or FORMAT number appears in a READ statement. 

RL Variable appears in a REAL statement. 

WR Variable or FORMAT number appears in a WRITE statement. 


29 


INOEX 


AAMAANANAAAAAAANANNAANAAAANAAAAAANAANAANDAAADNDANDANANAANAANAAAAANAANADNDIANDANAANAANANAAANAAADA 


a200 


PROGRAM -~ RAPLOT II 


PROGRAM -- RAPLOT II 
THE PURPOSE ‘OF THIS PROGRAM IS TO REDUCE THE RADIOACTIVITY SURVEY 
DATA FROM THE RIST PROJECT AND PLOT THE SURVEY ON A BENSON=LEHNER 
INCREMENTAL PLOTTERe THE FOLLOWING PLOTS ARE THE OUTPUT?® 
TRACKLINE FOLLOWED BY SURVEY VESSEL 
PLOT OF UNCORRECTED RADIATION VALUES (BACKGROUND SURVEY) 
SYMBOL PLOT OF RADIATION VALUES CORRECTED FOR BACKGROUND AND 
DECAY SINCE INJECTION TIMEe 
PROGRAMMER: PHILIP Ae TURNER 
GEOLOGY BRANCH 
US ARMY COASTAL ENGINEERING RESEARCH CENTER 
5201 LITTLE FALLS ROAD 
WASHINGTONe De Co 20016 
COMPLETED IN JANUARY 1969 


FORMAT AND ENTRIES ON DATA CONTROL CARD 
COL 1-3 CABLE LENGTH IN FEET TO THE NEAREST FOOT. 
COL 4= 6 WATER DEPTH PLUS FREEBOARD TO THE NEAREST FOOT. 
COL 7= 9 DISTANCE FROM RADAR MAST TO CABLE STANCHION IN FEET 
TO THE NEAREST FOOT. 
COL 10-16 ESTIMATED RACKGROUND COUNT RATE IN COUNTS/SEC FOR 
CHANNEL le 
COL 17-23 ESTIMATED BACKGROUND COUNT RATE IN COUNTS/SEC FOR 
CHANNEL 2 5 
COL 25-28 TIME OF INJECTION IN HOURS AND MINUTES. 
COL 29-31 THE NUMBER OF DAYS SINCE THE INJECTION 
COL 32-36 THE HALF LIFE OF THE RADIOISOTOPE IN DAYS. 
THE DECIMAL POINT MUST BE PUNCHED IN. 
COL 37-40 THE TIME WHEN THE SURVEY WAS STARTEDe IN HOURS AND 
MINUTESe : 
COL 41-42 TIME (IN SECONDS) BETWEEN SUCCESSIVE FIXES 
COL 44-50 BEACON1 NORTH COORD/ LAMBERT COORDINATES OF RADAR 
COL 52-58 BEACON1 EAST COORD/ BEACONS TO THE NEAREST FOOT. 
COL 60-66 BEACON2 NORTH COORD/ BEACONI IS ALWAYS UPCOAST. 
COL 68-74 BEACON2 EAST COORD/ 
COL 75-77 THE NUMBER OF LINES OF DATA TO BE SKIPPED AT THE 
BEGINNING OF A DATA SET IN ORDER TO AVOID 
READING IN SOME BAD DATA. 
FORMAT AND ENTRIES ON PLOT CONTROL CARD 
COL 1-3 PLOT OPTION CONTROL. TO USEe PUNCH THE NUMERAL 1 
IN THE COLUMN INDICATED. 
1 = PLOT TRACKLINE FOLLOWED BY SURVEY VESSEL. 


2 = PLOT UNCORRECTED RADIATION VALUESe 

3 = PLOT RADIATION VALUES CORRECTED FOR BACKGROUND 
AND DECAY SINCE TIME ZERO> 

4 = UNUSED. LEAVE BLANK. 


COL 5-14 MAP SCALE EXPRESSED IN FEET PER INCH 
COL 16-17 OPTION 1 / USE WHEN SPOTTING DATA FOR EACH PLOT 
COL 16-19 5) 2 7 OPTIONe USER CAN SPECIFY THAT EVERY NTH 
COL 20-21 J 4 / POINT BE PLOTTED. IF LEFT BLANKe THE 
PROGRAM ASSUMES EVERY POINT IS TO BE 
PLOTTED. 
COL 23-32 INTERVALS ON THE COORDINATE GRID AT WHICH TICK 
MARKS WITH THE LAMBERT COORDINATES WILL BE POSTEDe 
IF FIELO IS LEFT BLANKe PROGRAM WILL ASSUME THAT 
NO TICK MARKS ARE TO BE PLOTTED AND POSTED. 
COL 34-43 BEACON 1 / INJECTION SITE. DISTANCE IN METERS TO 
COL 45 54 BEACON 2 / THE NAMED BEACONS. IF FIELDS ARE LEFT 
BLANKe SITE IS NOT PLOTTED. 
COL 56-61 DAYeMONTH AND YEAR THE SAND WAS INJECTED 
FORMAT AND ENTRIES ON PLOT IDENTIFICATION CARD 
COL 1-78 THIS FIELD WILL BE PLOTTED ON THE LOWER MARGIN OF 
THE MAPo 
COL 80 PUNCH *°T*® HERE ON THE LAST DATA SET. 

COMMON NOPT(4) eSCALE*NPLT (3) pGRIDeSITENe SI TEEs INDATE oLEGEND (13) o BE 
LACLNe BEACIE e BEAC2Ne BEAC2E ep LINE BGCR (2) o SIGMA (2) eNBAR®EBAR eo NENTRY ON 
2MAXoNMINe EMAXe EMIN 

REAL NORTH(2200) »NCORD(2200) »NBARoNMAX oNMIN 

DIMENSION NMBR(2200)* TIME(2200)9 Dl2e2200)» EAST(2200)» ECORD(220 
10)+ RAD(202200)* CCR(202200) L(4)o IERR(2200)* BKG(2)° FATH(2200) 
2° ISIGN(3) 

INTEGER SENTNL(7) »WORD1°WORD2¢BLANK? AS 

LOGICAL JOBEND 


NTAPE IS THE NUMBER FOR THE INPUT TAPE CONTAINING THE SURVEY DATA 
ITAPE IS THE UNIT ON WHICH THE PROCESSED DATA IS WRITTEN OUT 


DATA WORD1»WORD2eNTAPEe ITAPEeCORLeCOR2¢BLANKe AS/GHEND OF 06H DATA o 
17°8020902e¢806H 0 1H*/ 
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PROGRAM == RAPLOT II 


READ .IN DATA CONTROL CARD 


IF (L(1)eNEe=1) GO TO 20 

GO TO 10 

IF (L(1)eLTe-1) GO TO 770 

READ (5¢780) CABLE* DEPTH» BOAT» BKGe ZHRe 2MINeDAYSeHLIFEo SET IME? RMING 
1SEC e BEACINe BEACLE e BEAC2Ne BEAC2E eo ISKIP 


READ IN PLOT CONTROL PARAMETERS 


READ (5°790) NOPT*SCALEe (NPLT(I) 0 T=193) eGRID»OUMP1 oDUMP2> INDATE 
DO 40 T=10e3 

IF (NPLT(I)) 30030040 

NPLT(T)=1 

CONTINUE 


READ IN PLOT LEGEND 


READ (50800) LEGENDe JOBEND 
WRITE (60810) LEGENDe JOBEND 
WRITE (60820) BEACINe BEACIE? BEAC2Ne BEAC2ZE 


COMPUTE PROGRAM PARAMETERS FROM DATA CONTROL CARD ENTRIES 


SQDSTB=(BEAC2N=BEACIN) ##2+ (BEACZE“BEACIE) ##2 
DISTa=SERT (S@DSTB) 

WRITE (6°830) SQDSTBrOISTB 
SINE=(BEAC2N=BEACIN) /DISTB 
COSINE=(BEAC2E=BEAC1E) /DISTB 

WRITE (6840) SINE? COSINE 

WRITE (6°850) ZHReZMINeSETIME + RMIN*’ SEC e DAYS 
ZHR=ZHR+ZMIN/60 6 

SETIME=SETIME+RMIN/606 
DELAY=SETIME+DAYS*24 .—-ZHR 

IF (HLIFE*GT0.0) DECAY=ALOG(20)/(HLIFE*24 0) 
WRITE (60860) HLIFE» DECAY» DELAY 

WRITE (69870) CABLE» DEPTH) BOAT 
CABLE=BOAT+SQRT (CABLE**2=DEPTH*#2) 

WRITE (6880) CABLE . 

WRITE (6890) (NOPT(T)eI=103) 


COMPUTE COORDINATES OF THE INJECTION SITE FROM THE DISTANCES FROM 
THE BEACONS 


IF (DUMP1) 80°80°50 
DUMP1=(DUMP1+COR1) #3228083 
DUMP2=(DUMP2+COR2) *3.28083 
O0X1=(SQDSTB+DUMP1*DUMP1-DUMP2*DUMP2) / (DISTB#2.0) 
OY 1=DUMP1*0UMP1=DX1*DXxX1 

IF (DY1) 60°60070 

SITEN=-999999. 

WRITE (60900) 

GO TO 80 

DY1==-SQRT(DY1) 
SITEE=DX1*COSINE-DY1*SINE+BEACIE 
SITEN=DX1*SINE+DY1*COSINE+BEACIN 
WRITE (60910) SITENeSITEE 


TEST FILE SENTINEL TO BE SURE THAT THE CORRECT DATA FILE IS BEING 
READ IN 


READ (NTAPE) SENTNL 

IF (SENTNL(1) eEQeWORDI ©ANDeSENTNL (2) eEQoWORD2) GO TO 730 

IF (SENTNL(1) eEQ@eLEGEND(12) eANDeSENTNL (2) EQ. LEGEND(13)) GO TO 90 
CALL NTRAN (NTAPE*601) 

GO TO 80 

WRITE (60920) SENTNL 


REAO IN THE DATA FILE FROM ONE RIST SURVEY 


TF (L(2)-NEo-1) GO TO 110 

6O TO 100 

IF (L(2)-LTo-1) GO TO 770 

CALL NTRAN (ITAPEo Le Ll 3SeLEGENDeL(1)) 

READ (NTAPE) LINEe (NMBR(N) eo TIME(N) eD( LON) oD620N) PRAD(CL9N) PRAD(20N) 
LoFATH(N) oN=10LINE) 

CALL NTRAN (NTAPE* 801) 

IERR(1)=0 

DO 130 N=2eLINE 
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PROGRAM == RAPLOT ITI 
TERR(N)=0 


CHECK TO SEE THAT LINE NUMBERS AND TIMES OF FIXES ARE IN MONOTONIC 
ASCENDING SEQUENCE. 


IF (NMBR(N) eLEoNMGR(N@1)) NMBR(N)=NMBR(NOL) +1 
IF (TIME(N)-TIME(N=1)) 12001200130 

TIME (N) =TIME(N=1) +SEC 

CONTINUE 

IF (ISKIPeLE«0) GO TO 160 


SKIP LEADING CARD IMAGES THAT CONTAIN BAD DATA. 


NSTART=ISKIP+1 

DO 150 N=NSTARTeLINE 
NMBR(N=ISKIP)=NMBR(N) 
TIME (N-ISKIP)=TIME(N) 
DO 140 I[=192 
D(IeN-ISKIP)=D(IeN) 
RAD(IeN-ISKIP)=RAD(T oN) 
FATH(N=ISKIP)=FATH(N) 
LINE=LINE-ISKIP 
MSTOP=LINE=1 


CHECK DISTANCES TO RADAR BEACONS FOR ERRORSe IF DISTANCE/TIME 
FOR SUCCESSIVE BEACON RANGES INDICATE A SHIP SPEED eGTe 6 KNOTS 
(3-09 METERS/SEC) RANGE IS IN ERRORe 


DO 270 I=le2 

DO 200 M=l1eMSTOP 

IF (D(IeM)) 2000200017C 
NSTART=M+1 

DO 190 N=NSTARTeLINE 

uF (ABS (O(ToN)—D( I 0M) DoC TIME (NI -TIME(M) )%3208865) 20072007180 
D(ieNd=—1e 
IERR(NJ=IERR(N) +1 
CONTINUE 

CONTINUE 

DO 260 M=1»MSTOP 

IF (D(IeM)) 21002100260 
NSTART=M 


CORRECT ERRONEOUS BEACON RANGES BY LINEAR INTERPOLATION (ON TIME) 
BETWEEN NON-ERRONEOUS RANGESe 


DO 230 N=NSTARToLINE 

IF (D(IeN)) 23002300220 

NSTOP=N 

GO TO 240 

CONT INUE 

OTIME=TIME (NSTOP) -TIME(NSTART=1) 

DD1=0( I eNSTOP)=D(I eNSTART=1) 
=NSTART 

D(IeN)=00I eNSTART=1) #DD01*% (TIME (N) =TIME(NSTART=1) ) /OTIME 
N=N+1 

IF (N=NSTOP) 25002607260 

CONT INUE 

CONTINUE 


COMPUTE POSITION OF SHIP FROM DISTANCES FROM THE TWO BEACONS 


LAG=0 
DO 300 N=leLINE 


MAKE CONSTANT CORRECTION FOR CUBIC AUTOTAPE INTERROGATOR 
AND CONVERT TO FEET 


OFT1=(D(19N) #COR1) #3228083 
DFT2=(D(20N)+COR2) #3228083 
DX1=(SQDSTB+DFT1*DFT1-DFT2*DFT2) /(2.*DISTB) 
DY1=OFT1I*OFT1I-DX12Dx1 


CHECK FOR IMAGINARY ROOT 


IF (DY1) 280°280°290 
TERR(N) =IERR(N) +4 
LAG=LAG+1 

GO TO 300 
OY1=-SORT(DY1) 
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adults COORDINATES AND TRANSLATE INTO CALIFORNIA LAMBERT COORDINAT 


EAST (N)=DX1*COSINE=DY1*SINE+BEACIE 
NORTH(N) =DX1*SINE+DY1*COSINE+BEACIN 
CONTINUE 

IF (LAG-EQ-0) GO TO 360 

ASSIGN 360 TO KEY 


ELIMINATE DATA SETS FOR WHICH AN 
IMAGINARY FIX WAS OBTAEFNED 


LAG=0 

DO 350 N=leLINE 

IF (IERR(N)=4) 330°3200320 
LAG=LAG+1 

GO TO 350 

NMBR (N=LAG) =NMBR (N) 
TIME (N=-LAG)=TIME(N) 
00 340 I=102 
O(TeN-LAG)=D(TeN) 

RAD (Ie N-LAG)=RAD(I oN) 
EAST (N=LAG) =EAST(N) 
NORTH(N=LAG)=NORTH(N) 
TERR (N=LAG)=IERR(N) 
CONTINUE 
LINE=LINE-LAG 

GO TO KEY? (3600430) 


CHECK NORTH AND EAST COORDINATES FOR 
EXTREME VALUES BY CHAUVENET*S CRITERION 


NBAR=-1. 

CALL STDEV (NORTHeLINE*?NBARe SONRTH) 
EBAR=-1. 

CALL STDEV (EASTeLINE*EBARe SDEAST) 
ALPHA=1e-1le/FLOAT (2*LINE) 
CHVR=TINORM( ALPHA’ $365) 

GO TO 370 

CHVR=5.0 

WRITE (60930) ALPHA 
GATE1L=EBAR-CHVR*SDEAST 
GATE2=EBAR+CHVR*SDEAST 
GATENIL=NBAR=CHVR*SDNRTH 
GATEN2=NBAR*+CHVR*SDNRTH 

LAG=u 

00 420 N=leLINE 

IF (EAST(N)-GATEL) 41073800380 
IF (EAST(N)=GATE2) 39003900410 
IF (NORTH(N)=GATEN1) 41004000400 
IF (NORTH(N)=GATEN2) 42004200410 
IERR(N)=IERR(N) +4 

LAG=LAGt1 

CONT INUE 

ASSIGN 430 TO KEY 


ELIMINATE ANY DATA SETS THAT HAVE AN EXTREME 
VALUES OF THE NORTH OR EAST COORDINATES 


IF (LAGeGTe0) GO TO 310 
CALL SUBROUTINE FOR PLOTTING THE TRACK OF THE SURVEY VESSEL. 


NENTRY=1 
CALL NTRAN (ITAPEeLeLINErNMBReL(2)) 
IF (NOPT(1)-EG@e1) CALL BENLH2 (NORTHe EAST oNMBR) 


APPLY A CORRECTION TO ALLOW FOR THE DISTANCE THE DETECTOR IS TOWED 
ASTERN OF THE SURVEY SHIPo 


IF (L(3)-eNEe-1) GO TO 450 

GO TO 440 

IF (L(3)eLTe-1) GO TO 770 
DNO=NORTH(1)=(NORTH(2)=-NORTH(1)) 
DEO=EAST(1)-(EAST(2)-EAST(1)) 

DENOM=SQRT ( (NORTH(1) -ONO) ##2+ (EAST (1) —-DEO) ¥#2) 
NCORD(1)=NORTH(1) -CABLE* (NORTH(1)=-DNO) /DENOM 
ECORD(1)=EAST(1)-CABLE* (EAST (1)-DEO) /DENOM 

DO 480 N=2eLINE 
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186 
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211 
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530 
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PROGRAM == RAPLOT II 
DENOM=SGRT ( (NORTH(N) =“NCORD(N=1) ) ®#2+ (EAST (N) -ECORD(N=1) ) ##2) 


THE CORRECTION FOR THE DISTANCE BETWEEN VESSEL AND THE DETECTOR 
IS EQUAL TO "CABLE*® UNLESS THE VESSEL IS LESS THAN *CABLE® FEET 


AWAY FROM THE LAST COMPUTED POSITION OF THE DETECTOR VEHICLE. IN 


THIS EVENTe THE NEW COMPUTED DETECTOR POSITION IS THE SAME AS 
THE LAST DETECTOR POSITION 


IF (DENOM=CABLE) 47004700460 

NCORD(N) =NORTH (NN) -CABLE*® (NORTH(N) “NCORD(N=1) ) /DENOM 
ECORD(N) =EAST(N) -CABLE* (EAST (N) -ECORD(N=1) ) DENOM 
GO TO 480 5 

NCORD(N)=NCORD(N=-1) 

ECORD(N) =ECORD(N-1) 

CONT INUE 

CALL NTRAN (ITAPEs Le LINE*NCORDeLAGe LoLINEsECORDeL(3)) 
IF (L(4)eNEe-1) GO TO 500 

GO TO 490 

IF (L(4).LTe-1) GO TO 770 


CONVERT RADIATION READINGS TO COUNTS PER SECOND 


IF (dKG(1)) 650° 650° 505 
CCR(201)=RAD(201)/SEC 
CCR(1°01)=RAD(101)/SEC 

00 510 N=2eLINE 

DO 510 I=1e2 
CCR(IeN)=RAD(I oN) /SEC 


COMPUTE THE MEAN AND STANDARD DEVIATION OF THE BACKGROUND COUNT 
RATE FROM THE RADIATION DATA THAT LIES WITHIN THE LIMITS OF THE 
ESTIMATED BACKGROUND COUNT RATE SET BY CHAUVENET®S CRITERION® 


IF (CCR(101) eGTe20*BKG(1)) CCR(101)=BKG(1) 
IF (CCR(201) eGTe20e#BKG(2)) CCR(201)=BKG(2) 
ALPHA=1 o-1 o/FLOAT (2*#LINE) 
CHVR=TINORM(ALPHA9 $515) 

GO TO S20 

CHVR=S.e 

WRITE (60930) ALPHA 

DO 590 I=le2 

BGCR(1I)=0-0 

SIGMA(T)=0.0 

NCOUNT=0 
GATE1=8KG(1)-CHVR*SQRT (BKG(1) /SEC) 
GATE2=8KG(1) +CHVR*SQRT (BKG(1)/SEC) 

DO 550 N=leLINE 

IF (CCR(IeN)=GATEL) 55005500530 

IF (CCR(IeN)=GATE2) 54005507550 
NCOUNT=NCOUNT+1 
BGCR(T)=BGCR(T)+CCR(ION) 

CONTINUE 

BGCR(I)=BGCR(I) /FLOAT (NCOUNT) 

DO 580 N=leLINE 

IF (CCR(I¢N)-GATE1) 58005800560 

IF (CCR(IeN)=GATE2) 57005807580 
SIGMA(I)=SIGMA(1)+(CCR(ILoN)-BGCR(I) ) #%2 
CONTINUE 
SIGMA(1)=SQRT(SIGMA( IS /FLOAT(NCOUNT) ) 
WRITE (60940) BKGeBGCReSIGMA 


CALL THE SUBROUTINE FOR PLOTTING UNCORRECTED RADIATION VALUES 
NENTRY=2 


IF (NOPT(2)cEQ@e1) CALL RDPLT2 (NCORDeECORDeCCR) 
IF (NOPT(3) oNEoi) GO TO 650 


CORRECT RADIATION VALUES FOR BACKGROUND COUNT RATE AND TIME=DECAY 


SUM=0.0 

SUM2=0.U 

NBAR=0-0 

EBAR=0.0 

00 620 N=1eLINE 

DO 610 I=1r2 

CCR(I»eN)=CCR(I oN) -BGCR(1) 

IF (CCR(I*ON) — 30*SIGMA(I)) 61096100600 

CCR(T oN) =(CCROI ION) =30*SIGMA(T) ) EXP (DECAY * (DELAY+TIME(N) /36000)) 


1 + 3.*SIGMA(I) 
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PROGRAM == RAPLOT II 


CONT [NUE 

IF (CCR(10N) eLE¢000) 60 To 620 
SUM=SUM4CCR(10N) 
SUM2=SUM2+CCR(20N) 


COMPUTE WEIGHTED MEAN AND STDe DEVe OF ACTIVITY LOCATION 


EBAR=E8aAR+(ECORD(N) -ECORD(1))*CCR(ieN) 

NB AR=NBAR+(NCORD(N) =NCORD(1)) *CCR(1L0N) 
CONTINUE 

NBAR=NCORD (1) +NBAR/SUM 
ESAR=ECORD(1)+EBAR/SUM 

SONRTH=0-0 

SDEAST=0.0 

00 640 N=1eLINE 

IF (CCR(1eN)) 64026400630 
SONRTH=SDNRTH+ (NEORD(N) -NBAR) ##2*%CCR(10eN) 
SDEAST=SDEAST+(ECORD(N) -EBAR) #*#2*CCR(10N) 
CONTINUE 

SDNRTH=SaRT (SDNRTHZSUM) 

SDEAST=SQRT (SDEAST/SUM) 

WRITE (60950) SUMeSUM2 

WRITE (6°960) NBAR*EBARe SONRTHe SDEAST 


COMPUTE AND PRINT 95 PCe CONFIDENCE LIMITS OF MEAN RADIATION 
LOCATION. 


RTSUM=SuRT (SUM/BGCR(1)) 

CF IDN=1-96*SDNRTH/RTSUM 

CF IDE=1.96*SDEAST/RTSUM 

WRITE (60970) CFION*eCFIDE 

DO 660 I[=ler2 

CCR(IeLINE+1)=BGCR(T) 
CCR(IeLINE+2)=SIGMA(T) 
NWRD=(LINE+2) #2 

CALL NTRAN (ITAPE*leNWRDeCCReL (4) ) 
NENTRY=3 

IF (NOPT(3)-E@e1) CALL RDPLT2 (NCORDeECORDe CCR) 
WRITE (60980) NMAXeEMAX 

WRITE (60990) NMIN*EMIN 


dRITe OUT THE NUMBERe COORDINATES AND ACTIVITY OF EACH DATA POINT 


XOUNT=50 

DO 720 nN=leLINE 

DO 690 J=1e3 

IF (FLD(3e-Jele TERR(N))) 67006700680 
ISEGN( J) =BLANK 

GO TO 690 

ISIGN(J)=AS 

CONT INUE 

IF (KOUNT=50) 71007000700 
WRITE _(6°1000) LEGEND 
KOUNT=0 


WRITE (6¢1010) NMBR(N) o TIME(N) oD(19N) oTSIGN(2) 0D (20N) oISIGN(2) 0 (RA 


1D(IeN) e L=102) pNORTH(N) PEAST(N) 0 ISIGN(3) eNCORD(N) eECORDI(N) o (CCR(TON 
2)eT=1e2) 


KOUNT=KOUNT+1 
WRITE (601020) 

IF (.NOTeJOBEND) GO TO 10 
CONTINUE 


ENO-FILE PLOT TAPE 


END FILE 9 

DO 760 I[=104 

IF (LCI) .NEo-1) GO TO 750 
GO TO 740 

IF (L(I)LTo-1) GO TO 770 
CONTINUE 

CALL NTRAN (ITAPE?9) 

STOP 


FORMAT (3F 300 2F7cO01X02F2c0eFSceOeF Sole SF2c004(1X0F 70) 013) 

FORMAT (4110F10-001Xe31201X0F100002(1X0F 10.0) 01 Xe AG) 

FORMAT (13A6eL2) 

FORMAT (10X013A6010XeL2) 

FORMAT (/5XeSHBEACON 1LeF10¢Oe1HNeF10eO*e1HEeSXeS8HBEACON 2°F10-00°1HN 


1eF10¢001HE//) 
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2c 

30 
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50 

60 

70 

80 

90 

log 
110 
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220 
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280 
290 
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316 
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365 
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FORMAT (5Xe21HSQUARE DIST BETWEEN =0616e8e¢5XeiBHDISTANCE BETWEEN = 
1°F10-0) 

FORMAT (5Xe6HSINE =0£160¢805XeB8HCOSINE =0E1608) 

FORMAT (5Xe1l6HINJECTION TIME =02F3-0010Xe12HCLOCK SET ATe2F3.005Xe 
AZ1HDIGITIZING INTERVAL =eF 300° THSECONDS/SXe30HDAYS ELAPSED SINCE I 
2NJECTION =0F 3-0) 

FORMAT (SXe22HHALF=LIFE OF ISOTOPE =eF7e2e4HDAYSeSXe14HDECAY FACTO 
IR =0£1668°SXeLSHTIME=DELAY FACTOR =eF7e2*eSHHOURS) 

FORMAT (10X°e14HCABLE LENGTH =eF4e0°010Xe18HMEAN WATER DEPTH =eF4.00 
110X°13HBOAT LENGTH =eF4.0) 

FORMAT (10X¢38HDISTANCE FROM RADAR MAST TO DETECTOR =eF6.10SHFEET. 
1) 

FORMAT (//40Xe15HPLOTS GENERATED/10Xe SHTRACKLINEe 150 10Xe 20HBACKGRO 
LUND RADIATIONe®I5010X*°19HCORRECTED RADIATION? IS) 

FORMAT (5Xe51HBEACON RANGES FOR DUMP SITE COMPUTE IMAGINARY ROOT.) 

FORMAT (10Xe37HLAMBERT COORDINATES OF INJECTION SITE*F10.001HNeF10 
1.0°1HE) 

FORMAT (1HO019X0 7A6) 

FORMAT (37H THERE WAS AN OVERFLOW WHEN ALPHA WASoF6e3eSXe2SHCHVR W 
AAS SET EQUAL TO 5-0) 

FORMAT (//20X¢S3HSUMMARY STATISTICS OF BACKGROUND RADIATION COUNT 
ARATE/30Xe13HRAD CHANNEL 105Xe13HRAD CHANNEL 2/10Xe21HEST. RBKGe COU 
2NTS/SECo oF 10000 8X0F1L0c00/10X°21HMEAN BKGe COUNTS/SECe oF10.008X0F 106 
BO/1OXP2LHSTO*e DEVe COUNTS/SECe 0 F100008X0F 1000) 

FORMAT (//20X°33HSUM OF CORRECTED RADIATION COUNTS/25X»0 19HRADIATIO 
IN CHANNEL 1°E16e8/25Xe19HRADIATION CHANNEL 2°E£16-8) 

FORMAT (//20Xe41HSUMMARY STATISTICS OF RADIATION LOCATION. /24Xe11H 
ANORTH COORD?10Xe1LOHEAST COORD/16Xe YHMEANe 5X0F10-0010X0F10-0/11X 09H 
2STDe DEVerSXeF10¢0010X0F 10-0) 

FORMAT (10XeLOHCONFIOENCE/7Xe13HLIMIT OF MEANeSXeF10-0010X0F10-0) 

FORMAT (//7Xe13HMAXIMUM COORDe5Xe0F 10000 1LHNe9XeF 10-00 1HE) 

FORMAT (//7Xe13HMINIMUM COORDeSXeF10¢0e1HNe9XeF10e001HE/S/) 

FORMAT (1H109X013A6//2Xe1OHLINE TIME0SXe1BHDISTANCE TO BEACONeSXoe 
LISHRADIATION COUNT?¢4Xe16HBOAT COORDINATES*4Xe16HBALL COORDINATES?3 
2Xe 1QHCORRECTED RADIATION/9Xe SHSECe 11X09 1H1 9 10X09 1H206XeSHRADI 6X 4HR 
BAD20 9X0 LHNe 9X0 LHEe 9X0 LHNe 9X0 LHE*e 7XeGHRADI» 3X0 9H RAD2) 


FORMAT (1XeTISeF6cOoFileleALoF1OcleAlo4F1O-OrAle2Fi0.O0F11.00F11.0) 
FORMAT (1H1) 
END 
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83 B4% 
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91 92 95" 
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99" 101 
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116 119 135 136% 
143% 
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EXP 
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GRID 
HLIFE 
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IERR 
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J 
JOBEND 
KEY 
KOUNT 
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LEGEND 
LINE 
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248WR 
243 


80 
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253 
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98= 
lal 
151 
178 
213 
240 


NBAR 
NCORD 


NCOUNT 
NENTRY 
NMAX 
NMBR 
NMIN 
NOPT 
NORTH 


NPLT 
NSTART 
NSTOP 
NTAPE 
NTRAN 
NWRD 
RAD 
RDPLT2 
RMIN 
RTSUM 
SCALE 
SDEAST 
SDNRTH 
SEC 
SENTNE 
SETIME 
SIGMA 
SINE 
SITEE 
SITEN 
S@DSTB 
SORT 


STDEV 
STOP 
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SuUM2 
TIME 
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241 242 243 263 265 273WR 
ico 2RL 136= 137AG 147 148 
248WR 
2Ru 169= 172 174= 177= 180AG 
259AG 273WR 
200= 206= 209 215 
1co 160= 217= 258= 
ico 2RL 26008R 
301 58RD 63 70= 125= 161AG 
1c0 2RL 261WR 
ico 11RD 34wR 162 218 219 
eRe 116= 131= 137AG 153 154 
172 174 273wR 
1co 11RD 13 14> 
68= 69 81= 62 90= 91 
93= 96 97 101 
6DA 48RD 51AG 58RD 59AG 
51 57 59 161 180 257 
256= 257AG 
301 58RD T4s 129= 185 186 
216 259 
10RD 25wWR 27 
249= 250 251 
1co 11RD 
139AG 145 146 239= 243= 246= 
137AG 147 148 238= 242= 245= 
10RD 25wWR 65 185 186 169 
GIN 48RD 49 50 S3wR 
10RD 25wWR 27= 28 
1C0 199= 213= 215= 216WR 227 
22= 24wR 4S 46 115 116 
1c0 45= 47WR 
1C0 41= 46= 47WR 
19= 20 21wR 38 108 
20 32 4y 114 168 172 
246 249 
137 139 
285 
220= 231= 236 | 237 245 246 
221= 232= 247WR 
301 58RD 64 65= 71= 83 
273WR 
141 193 
GIN 6DA 49 
GIN 6DA 49 
10RD 25wWR 26= 28 
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SUBROUTINE BENLH2 (NORTHEAST oe NMBR) 


SUBROUTINE BENLH2 (NORTH? EAST NMBR) 
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273WR 
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97 


273WR 


251 
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202 


255 


202 


249 
99 


236= 


236 


168 


98 


215 
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BoL 
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THIS SUBROUTINE GENERATES THE PLOT INSTRUCTIONS TO DRIVE A BENSON-B=-L 
LEHNER MODEL 305 DIGITAL PLOTTER. IF THE INSTRUCTIONS ARE BEING 8-L 


WRITTEN ON TAPEr USE OPTIONS Ne Ke AND E ON THE ASSIGN CARDe 


PROGRAMMER’ PHILIP Ae TURNER»e GEOLOGY BRe U S A CERC 


BeL 
BeL 
Be-L 


COMMON NOPT(4) eSCALE*NPLT(3) oGRID*DUMPNe OUMPE e INDATE s LEGEND (13) eBEB-L 
LACAN? BEAC1Es BEAC2Ne BEACZE o NPTSeBGCR (2) SIGMA (2) eNBAReEBAReNENTRY oNB-L 


2MAX eNMINeEMAXe EMIN 


B= 


DIMENSION EAST(2200)° RAD(2°2200) NMBR(2200)0 x(2200)° Y(2200)» RA-L 


AADUL (7) 
REAL NORTH(2200) eNBAReNMAXeNMINeNORGIN 
ITRNP=0 
GO TO 10 
ENTRY RDPLT2 (NORTH*eEAST*RAD) 
IF (NOPT(1)-EQ@e1) GO TO 120 
IF (NOPT(2) eEQ@e1LeANDeNENTRYeEQe3) GO TO 120 


BeL 
Bo-L 
BeL 
BeoL 
BeL 
BeL 
BeoL 


DEFINE THE LOGICAL UNIT ON WHICH PLOTTER COMMANDS ARE TO BE OUTPUTAR=L 


CARD PUNCH = 3 
MAGNETIC TAPES = 7°34 


CALL INPLOT (9) 


(ail 
BeL 
Bet 
BeL 
BaL 


DETERMINE THE MAXIMUM ANO MINIMUM VALUES FOR BOTH COORDINATES~ B= 


39 


ReL 


242 
242 


169 


99 


245 


228 


w~OrINoufuUNre 


ANNONh 


SUBROUTINE BENLH2 (NORTHEAST? NMBR) 


NMAX=NORTH(1) 
NMIN=NORTH(1) 


EMAX=EAST (1) 
EMINSEAST(1) 


DO 20 N=2eNPTS 

NMAX=AMAX1 (NMAXe NORTH(N) ) 
NM IN=AMINL (NMINeNORTH(N) ) 
EMAX=AMAX1 (EMAX* EAST(N) ) 
EMIN=AMIN1 (EMINeEASTI(N) ) 
CONTINUE 


CHECK SCALE TO BE SURE THAT NO FIXES WILL PLOT OFF THE MAPe 
CHANGE SCALE UNTIL ALL THE FIXES FIT ON THE PLOT. 


YSCALE=SCALE 


IF (NMAX=(NMIN+25¢*SCALE)) 60°60°30 
ITRNP=1 
IF (EMAX=(EMIN+25¢*SCALE)) 60060040 
IF ((EMAXEMIN) =(NMAX=NMIN)) 60050050 
ITRNP=0 
IF (ITRNPcEQ@e0) GO TO 90 
ASSIGN 90 TO KEY 


CELL= 


EMAX 


EMAX=NMAX 
NMAX=-EMIN 
EMIN=NMIN 
NMIN==CELL 
DO 80 N=leNPTS 
CELL=NORTHI(N) 
NORTH(N)=-EAST(N) 


EAST(N) =CELL 


GO TO KEYe (900130) 
IF (NMAX=(NMIN+25¢*SCALE)) 11001100100 
SCALE=SCALE+YSCALE 
GO TO 90 


BL 
BeL 
Bei 
Bel 
BeaL 
RL 
B-L 
ReaL 
BeL 
B=" 
B=G 
B= 
Bae 
Be-L 
Bel 
BeL 
BeoL 
BeL 
BeL 
Bol 
Bel 
BeL 
BeL 
BeL 
BeoL 
B= 
Be=L 
BeL 
BaL 
B= 
B-L 
BaL 
BeL 
BeoL 
BeoL 
(AL 


TAKE THE ORIGIN OF THE PLOT AS THE MINIMUM COORDINATES REMAINDEREDB-L 
BY THE PLOT SCALEo 


NORGIN=SCALE*AINT (NMIN/SCALE) 
YORGIN=NORGIN/SCALE 

EORGIN=SCALE*AINT (EMEN/SCALE) 
XORGIN=EORGIN/SCALE 

IF (NENTRYeEQ@e3eANDeNOPT(2) eEQe1) GO TO -150 


OIVIDE THE COORDINATES (IN FEET) BY THE SCALE (IN FEET). 


IF (NENTRYcEQ@oleOReITRNPcEQe0) GO TO 130 
ASSIGN 130 TO KEY 

GO TO 70 
DO 140 N=1eNPTS 

X (NW) SEAST(N) SCALE 
¥ (N) =NORTH(N) /SCALE 


WRITE IN THE LEGENDe SCALE* TYPE OF PLOT AN INJECTION DATE AT THE 


BOTTOM LODGE 


CALL 
CALL 
CALL 
CALL 


LETTER 
LETTER 
NUMBER 
LETTER 


OF THE PLOT. 


(78080001¢00—-1e0*LEGEND) 
(8080005e00-2e2°8H1 INCH =) 
(SCALE0 40108000 Te20~20e2) 
(4080009e80-2.2°4HFEET) 


GO TO (16001700180)° NENTRY 


CALL 


LETTER 


GO TO 190 


CALL 
CALL 


LETTER 
LETTER 


GO TO 190 


CALL 
CALL 
CALL 
CALL 
CALL 
CALL 


LETTER 
LETTER 
LETTER 
LETTER 


(25080001200—202025HPLOT OF SURVEY TRACK LINE) 


(29060001200-2e0e29HPLOT OF UNCORRECTED RADIATION) 
(29960001200-2e6*29HROUNDED TO NEAREST 100 COUNTS) 


(28060001200-2e0e28HPLOT OF CORRECTED RADIATION/) 
(27060001200—2e6727HBACKGROUND RADIATION RATIOS) 
(140600020¢0-2e0014HINJECTION DATE) 
(60600021080-2-50 INDATE) 


PLOT (00000003) 
PLOT (XORGINe YORGING1) 


40 


B-L 
BeL 
BeaL 
BeL 
Bel 
ReL 
BeL 
BeL 
BeL 
B= 
B-L 
BeL 
BeL 
BeL 
BeL 
B=6 
B=— 
Bets 
Rel 
BeL 
BeL 
BeL 
BeL 
Bab 
BeL 
BL 
BeL 
Bel 
BeL 
Balk 
BeL 
Bab 
BeL 
BeL 
BeL 
Bel 
B-L 


SUBROUTINE BENLH2 (NORTHeEASTeNMBR) 


c PLOT IN TIC MARKSe WITH COORDINATES» AT THE LEFT AND LOWER BORDERSA-L 100 
Cc OF THE PLOT. B-L 101 
Cc RB-i 102 
69 CALL SYMBOL (XORGINe YORGINe 20066) B-L 103 
70 IF (GRIDeLE»000) GO TO 230 B-L 104 
71 TICN=GRID* (10 +AINT(NORGIN/GRID) ) B-L 105 
72 TICE=GR10*(1.0+AINT(EORGIN/GRID) ) B-L 106 
73 YTIC=TICN/SCALE R-L 107 
74 XSHIFT=XORGIN=1.28 A-L 108 
75 200 CALL REPSYM (XORGIN®YTIC) B-L 109 
76 CALL NUMBER (TICN*8e004eO0eXSHIFTYTIC) B-L 110 
77 YTIC=YTIC+GRID/SCALE Bol 111 
78 TICN=TICN+GRID B-L 112 
79 IF (YTIC=(YORGINt+25-)) 20002000210 B-L 113 
80 210 =XTIC=TICE/SCALE Bei 114 
b1 YSHIFT=YORGIN-e5 B-L 115 
82 220 CALL REPSYM (XTIC*YORGIN) B-L 116 
83 CALL NUMBER (TICE*800°4e00XTICe YSHIFT) B-L 117 
84 XTIC=XTIC+GRID/SCALE B-L 118 
85 TICE=TICE+GRID 5 B-L 119 
86 IF (XTIC=(XORGIN+25¢)) 22002209230 B-L 120 
87 230 CONTINUE B-L 121 
Cc B-L 122 
Cc PLOT THE INJECTION SITEe B-L 123 
Cc R-L 124 
88 IF (OUMPNeGTeNMAXeOReDUMPNeLTeNMIN) GO TO 240 B-L 125 
89 YD=DUMPN/SCALE B-L 126 
90 XD=DUMPE/SCALE B-L 127 
91 CALL SYMBOL (XDeYDe10070) B-L 128 
92 XD=XU+02 B-L 129 
93 CALL LETTER (49490°XD*YDe4HDUMP) B-L 130 
Cc Bol 131 
Cc PLOT THE MEAN LOCATION OF THE RADIATION DISTRIBUTION. B-L 132 
c B-L 133 
94 240 IF (NENTRY=3) 260°250°260 RB-L 134 
95 250 IF CITRWPeEQ@e1) GO TO 260 R-L 135 
96 YB=NBAR/SCALE A-t 136 
97 XB=EBAR/SCALE B-L. 137 
98 CALL SYMBOL (XBeYBe10071) B-L 138 
99 XB=XB+02 B-L 139 
100 CALL LETTER (6¢4°0°XBr Y8e6HRADBAR? B-L 140 
101 260 NN=NPLT(NENTRY) BoL 141 
102 NSTART=14NN B-L 142 
103 IF (NENTRY=2) 27002900310 R-L 143 
Cc B-L 144 
c PLOT THE TRACK LINE FOLLOWED BY THE SURVEY VESSELe B-L 145 
C B-L 146 
104 270 =©CALL LINE (Xe YoNPTSe306505) Rel 147 
105 DO 280 N=1eNPTSeNN B-L 148 
106 Z=FLOAT (NMBR(N) ) Bel 149 
107 CALL NUMBER (Z93000200°X(N) eY(N)) B-L 150 
108 280 CONTINUE B-L 151 
109 GO To 440 ReL 152 
C R-L 153 
c PLOT THE BACKGROUND RADIATION AT EVERY NNTH POINT Bal 18 
c B= 
110 290 CALL SYMBOL (X(1)0Y(1) 03067) B-L 156 
111 Z=RAD(101)/100. R-L 157 
112 CALL NUMBER (Z030002000X(1) oY(1)) R-L 158 
113 DO 300 I=NSTARTeNPTS B-L 159 
114 CALL REPSYM (X(1) 0¥(1)) B-L 160 
115 IF (HOD(IeNN) eNE©O) GO TO 300 Bel 161 
116 Z=RAD(10T) 71006 Bel 162 
117 CALL NUMBER (Z030002000X(1) 0 Y(I)) R-L 163 
118 300 CONTINUE B-L 164 
119 GO TO 440 B-L 165 
Cc B-L 166 
Cc PLOT CORRECTED RADIATION VALUES AT EVERY NNTH POINTe at ey 
Cc Re 6 
120 310 RADUL(1)=3.*SIGMA(1) A-L 169 
121 RAUUL(2)=1006 A-L 170 
122 00 320 L=3e7 B-L 171 
123. 320 RADUL(L)=2.*RADUL(L=1) Re-L 172 
124 LAST=0 R-L 173 
125 DO 420 I=1leNPTS R-L 174 
Cc A-L 175 


4 


i2o 
127 
128 
129 
130 
131 
132 
135 
134 
155 
1s5 
137 
138 
139 
140 
141 
142 
445 
44 
145 
146 
147 
148 
149 
150 
151 
152 
153 
154 
155 
156 
157 
158 
159 
160 
161 
162 
163 
164 
165 


166 
167 
168 
169 
170 
171 
172 
173 
174 
175 
176 
177 
178 
179 
180 
181 


SYMBOL. 


10 
20 
30 
40 
50 
60 
70 
80 
90 
100 
110 


NON0a0 


330 


340 


350 
360 
370 
380 
390 
400 
410 


420 


430 


450 
460 


SUBROUTINE BENLH2 (NORTH*eEASToNMBR) 


DATA POINTS HAVING CORRECTED RADIATION COUNTS MORE THAN 3 STANDARDA-L 


DEVIATIONS BELOW MEAN BACKGROUND COUNT RATE ARE NOT PLOTTED. 


SUCHR=-L 


LOw READINGS MAY INDICATE THAT THE DETECTOR WAS OVERTURNED OR WAS A-L 


"FLYLNG® 


IF (RAD( Le T)+RADUL(1)) 42004200330 
JO 350 L=107 

IF (RAD(1eT)=RADUL(L)) 340°3400350 
NSYM=644L 

Ir=1 

IF (ySYiteGTeo5) TH=2 

IF (NSYm%eGTeb9) ITH=4 

GO TO 300 

CONTINUE 

NSYM=72 

TH=4, 

IF (NSYM=LAST) 380¢3700386 

CALL REPSYM (xX(i) 0 V(I)) 

GO TO 390 

CALL SYMBOL (X(T) 0 Y(1) eIHe NSYM) 
LAST=NSYM 

iF (MOO(ITeNN)) 400740004206 

IF (RAD(1LeT)-BGCR(1)) 42004200410 
Z=RAD(1eT)/BGCR(1) 

CALL NUMBER (Ze30002rO0eX(I)eVC(I)) 
CONTINUE 2 
XEDGE=AMAX1 ( (EMAX/SCALE) +2. 0 XORGIN+27) 
YPT=YORGIN+1205 


CALL LETTER (2604°00XEDGEs YPTe26HCORRECTED COUNT RATE (C/S)) 


YPT=YPT--5 

XSHIFT=XEDGE+ eS 

XSHF T2=XSHIFT+200 

NSYM=64 

IH=1 

DO 430 L=1e7 

NSYM=NSYM*1 

IF (NSYMeGTe65) IH=2 

IF (NSYMeGT2e69) IH=4 

CALL SYMBOL (XEDGE* YPT» IHeNSYM) 


CALL LETTER (1603000XSHIFT*YPT*e16GHCCR ol Te 


CALL NUMBER (RADUL(L) 05¢1¢3°0°eXSHFT20YPT) 
YPT=YPT—04 

CALL SYMBOL (XEOGE* YPTe THe 72) 

CALL LETTER (16¢3000XSHIFTo YPT»16H 

CALL NUMBER (RADUL(7) 0501030 00XSHFT20YPT) 


MOVE PEN TO RIGHT BORDER OF PLOT IN PREPARATION FOR NEXT PLOT. 


XEDGE=AMAX1 (EMAX/SCALE49 0 op XORGING3S4 o ) 
CALL PLOT (XEDGE* YORGINe3) 
CALL PLOT (0e000c001) 

CALL PLOT (0-000¢00=3) 

IF CITRNP°EQ@-0) GO TO 460 
CELL=EMAX 

EMAX=-NMIN 

NMIN=EMIN 

EMIN=-NMAX 

NMAX=CELL 

DO 450 N=1leNPTS 
CELL=EAST(N) 
EAST(N)==NORTH(N) 

NORTH(N) =CELL 


RETURN 

END 

6 10* 

15 20* 

22 23* 

24 25% 

25 26% 

22 24 25 27% 
29% 49 

34 37* 

BY 28 38 39% 41 
39 4O* 

39 42x 


42 


GAL 
Real 
BeL 
Beat 
ReL 
BL 
B-L 
Bel 
B= 
BeL 
Bel 
BeaL 
Rak 
BeL 
Bel 
Bel 
Rel 
Bel 
Rel 
BeL 
BeL 
Bel 
Bel 
BeL 
BeL 
BeL 
BeL 
BeL 
Bel 
BeL 
Bel 
BeL 
BeL 
B-L 
BeL 
BeoL 
BeL 
BeL 
BeL 
BeL 
BeL 
B= 
BeL 
BeL 
BeL 
BeL 
Bol 
BeL 
BeL 
BeL 
BeoL 
Bel 
BoL 
BeL 
BeL 
BeL 
BeL 
BeL 
BoL 
BoL 
BeL 


176 
177 
178 
1719 
180 
181 
182 
183 
184 
185 
186 
187 
188 
169 
190 
191 
192 
193 
194 
195 
196 
197 
198 
199 
200 
201 
202 
203 
204 
205 
206 
207 
208 
209 
210 
211 
212 
213 
214 
215 
216 
217 
218 
219 
220 
221 
222 
223 
224 
225 
226 
227 
228 
229 
230 
231 
232 
233 
234 
235 
236 
237 
238 
239- 


SUBROUTINE BENLH2 (NORTH? EAST oNMBR) 


120 - 8 9 46x 

130 - 38 47 48 50% 

140 - 50 52x 

150 - 46 53% 

160 - 57 58* 

170 - 57 60% 

180 - 57 63% 

190 - 59 62 65% 

260 - 75* 79 

210 - 79 80* 

220 - 82% 86 

230 - 70 86 87% 

240 - 88 94x 

250 - 94 95% 

260 - 94 95 101 

270 = 103 104% 

280 - 105 108* 

290 = 103 110* 

300 = 113 115 118% 

310 - 103 120% 

320 - 122 123* 

330 = 126 127% 

340 - 128 129% 

350 - 127 128 134% 

360 = 133 137% 

370 = S72, 138% 

380 - 137 140% 

390 - 139 142% 

400 = 142 143% 

410 = 143 144% 

420 =) 125 126 142 143 146% 

430 = 155 162* 

440 = 109 119 166* 

450 - 176 179* 

460 = 170 180* 

AINT - 42 44 71 72 

AMAXL == 16 18 147 166 

AMINL) = at? 19 

BEAC1IE = 2co 

BEACIN = 2co 

BEAC2E = 2co 

BEAC2N = 2co 

BENLH2 = 1 

BGCR - 2co «143 144 

CELL - 29= 33 35= 37 171= 175 - 177= 179 

OUMPE = 2co 90 

DUMPN = 2co 88 89 

EAST - 1AG 3o1 7 13 14 18 19 36 37= 51 
177 178= 

EBAR - 2co 97 

EMAX - 2co 13= 18= 24 25 29 30= 147 166 171 
172= 

EnIN - 2co 14= 19= 24 25 31 32= a4 173 174= 

EORGIN = 44s 45 72 

FLOAT = 106 

GRID - 2co 70 71 72 77 78 64 85 

I oO 113 114AG 115 116 125 126 128 138AG 140AG 142 
‘143 144 

TH = 130= i31= 132= 136= 140A6 154= 157= 156= 159AG 163AG 

INDATE = 2co 

INPLOT = 10 

ITRNP = 5= 23= 26= 27 47 95 170 

KEY = 28= 38 4g= 

L = 122 123 127 128 129 155 161AG 

LAST = 1245 137 141= 

LEGEND = 2co 

LETTER = 53 54 56 58 60 61 63 64 65 66 

93 100 149 160 164 

LINE = 104 

MOD - 115 142 

N = 15 16 17 18 19 34 35 3e 37 50 

Sr 52 105 106 176 177 178 179 

NBAR - 2co QRL 96 

NENTRY = 2co 9 46 47 57 94 101 103 

NMAX - 2co GRE 1l= 16= 22 25 30 31= 39 88 
174 175= 

NMBR - 1AG 301 106 

NMIN - 2co 4RLE 12= 17= 22 25 32 33= 39 42 
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SUBROUTINE BENLH2 (NORTHEAST eNMBR) 


88 172 173= 
NN = 101= 10 105 115 142 
NOPT S 2co 8 Re a6 
RGIN = SRL 42s 
AGERE = 1AG SRL 7 11 12 16 17 35 36= 52 
178 179= 
T = 2co 101 
here = 2co 1s 34 50 109A6 105 113 125 176 
NSTART = \02= 113 
NSYM a 129= 131 132 135= 137 140AG 141 153= 156= 157 
158 159AG 
NUMBER = 55 76 83 107 112 117 145 161 165 
PLOT = 67 68 167 168 169 
RAD = 301 U 111 116 126 128 143 144 
RADUL o SDI 120= 121= 123= 126 128 161AG 165AG 
RDPLT2 = Uf 
REPSYM = 75 82 114 138 
RETURN = 180 
SCALE — 2co 21 22 24 39 40= 4e 43 ay 4S 
Sl 52 SS5SAG 73 77 80 64 89 90 96 
97 147 166 
SIGMA = 2co 120 
SYMBOL = 69 91 98 110 140 159 163 
TICE Oo 72= 80 83AG 85= 
TICN = 71l= 73 76AG 78= 
x = 301 Sl= 104AG 110AG 114AG 138AG 140A6 
XB cS 97= 9BAG 99= 
xn = 90= 91AG 92= 
XEDGE o 147= 151 159AG 163AG 166= 167AG 
XORGIN = 45= 68AG 69AG 74 7SAG 86 147 166 
XSHFT2 = 152= 
XSHIFT = 74= 151= 152 
xTICc o 80= 82AG 84= 86 
Y = 3DI1 52= 104AG 110AG 114AG 138AG 140AG 
YB ° 96= 9BAG 
Yo on 89= 91AG 
YORGIN <= 43= 68AG 69AG 79 61 82AG 148 167AG 
YPT o 146= 150= 159AG 162= 163AG 
YSCALE =- 21s 40 
YSHIFT = 61= 
yTic = 73= T5AG 77= 79 
2 oO 106= 107AG 111= 112AG 116= 117AG 1445 145AG 
1 SUBROUTINE STDEV(XoNeXNDeS) STOEV 
(eer ween nnn nner enero eee e nto e se ere tee ec eens eeesoresasesesososesesesese=STDFV 
Cc CALCULATES THE STANDARD DEVIATION OF A SEQUENCE OF DATA POINTS STDEV 
(eet ween nnn nnn eee etn ee cece tees ee setae teset ese sseesesesesesesesssosore=sSTOFV 
2 OIMENSION X(1) STDEV 
3 ENS=XND STDEV 
4 CALL AMEAN(XeNeXND) STDEV 
5 S=00 STOEV 
6 DOL L=10N STDEV 
U 1 S=S4(X(1)—-XND) #(X (1) =XND) STDEV 
8 TF(ENS eLTe Oo) GO TO 2 STDEV 
9 N=N-1 STDEV 
10 2 S = SQRT(S/FLOATI(N)) STDEV 
il RETURN STDEV 
12 END STDEV 
SYMBOL PRBSeBeSeSepeesagss REFERENCES BOSBSBRB BBS Ssessesss 
1 = 6 7* 
2 oO 8 10* 
AMEAN - 4 
ENS = 3= 8 
FLOAT = 10 
I 2 6 7 
N = 1AG 4AG 6 9= 10 
RETURN =- 11 
Ss = 1AG 5= 7= 10= 
SORT o 10 
STOEV = 1 
x = 1AG 201 GAG Uv 
XND = 1AG 3 4AG 7 
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ONOuUFUN 


SYMBOL 


1 
AMEAN 
FLOAT 
I 


N 
RETURN 
x 

XBAR 


1 
2 
3 


4 
5 
6 
v 
8 
9 

10 

ll 

12 

13 

14 


SYMBOL 


1 

A 

ALOG 
ALPHA 
B 

I 
OVERFL 
RETURN 
SQRT 
TINORM 
x 


SUBROUTINE AMEAN(XeNe XBAR) 


SUBROUTINE AMEAN (Xe Ne XBAR) AME AN 
(GOO ESS Fa a ee ee ew awn re ——AMEAN 
Cc CALCULATES THE ARITHMETIC MEAN OF A SEQUENCE OF DATA POINTS AMEAN 
(Sr ee ee een n ee tae sete toto seonres=—AMEAN 

DIMENSION X(1) AMEAN 

XBAR=06 AME AN 

001 I=1leN_ |. AMEAN 

1 XBAR=XBAR+t+X(1) AME AN 

XBAR = XBAR/FLOAT(N) AMEAN 

RETURN AMEAN 

END AMEAN 

SES8BS68 888 e2oe8 es 8 REFERENCES PESBSBSBSSsesesee ses 
4 5* 
1 
6 
4 5 
1AG 4 6 
7 
1AG 201 5 
1AG 3= 5= 6= 

FUNCTION TINORM(ALPHA?S) TINORM 

DIMENSION A(3) 0B(3) TINORM 

DATA(A(T) e T=103) 700103280 .802853020515517/0 (B(1) e1=103)/.0010308* TINORM 

1. TORE D0 BORE VEY TINORM 
(BESS SE Sa Sear arama. eee terete sete se sesece===T INORM 
Cc [APPROXIMATION TO INVERSE NORMAL DISTRIBUTION TINORM 
Cems e tes eee tease a= ee ee ete fetes erases esesesesese2es2=T INORM 

IF (sNOTs (ALPHAsGTeO+sANDeALPHAcLTole)) Go TO 1 TINORM 

X=ALPHA TINORM 

IF (XeGTe0e5) X=10e-X TINORM 

X=SORT(-20 *ALOG(X)) TINORM 

TINORM=X=—(A(3)4X#(A(2)4X9A(1)))70104X*(B(3)+X8(B(2)4xX* B(1)))) TINORM 

CALL OVERFL(I) TINORM 

IF(1eEQ@e1) RETURN 2 TINORM 

IF(ALPHAcLTee5) TINORM=-TINORM TINORM 

RETURN TINORM 

1 RETURN 2 TINORM 
END TINORM 
PSseseesgesa8s58 sag see REFERENCES SE88S8S88SeeseRB ese 88 
4 13* 
2DI 3DA 8 
7 
1AG 4y 5 11 
201 3DA 8 
3DA 9AG 10 
9 
10 12 13 
7 
1 = lis 
5= 6 7= 8 
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APPENDIX B 


LISTING AND INDEX OF RAPLOT III PROGRAM 
Ee SASS 


FORTRAN IV Listing of RAPLOT III and Subroutines 
TRACK and TINORM with an Index to all Statement 
Numbers, Variable Names, and Subroutine Calls 
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PROGRAM -- RAPLOT III MODIFIED FOR IBM-7094 WITH SC-4060 CRT 


PROGRAM == RAPLOT III MODIFIED FOR IBM=7094 WITH SC=4060 CRT 
THE PURPOSE OF THIS PROGRAM IS TO REOUCE THE RADIOACTIVITY SURVEY 
DATA FROM THE RIST PROJECT AND PLOT THE SURVEY ON A BENSON@LEHNER 
INCREMENTAL PLOTTERe THE FOLLOWING PLOTS ARE THE OUTPUT? 
TRACKLINE FOLLOWED BY SURVEY VESSEL 
PLOT OF UNCORRECTED RADIATION VALUES (BACKGROUND SURVEY) 
SYMBOL PLOT OF RADIATION VALUES CORRECTED FOR BACKGROUND AND 
DECAY SINCE INJECTION TIME. 
PROGRAMMERS PHILIP Ae TURNER 
GEOLOGY BRANCH 
US ARMY COASTAL ENGINEERING RESEARCH CENTER 
5201 LITTLE FALLS ROAD 
WASHINGTONe De Co 20016 
COMPLETED IN JANUARY 1969 


FORMAT AND ENTRIES ON DATA CONTROL CARD 
COL i= 3 CABLE LENGTH IN FEET TO THE NEAREST FOOTe 
COL 4= 6 WATER DEPTH PLUS FREEBOARD TO THE NEAREST FOOTe 
COL 7= 9 DISTANCE FROM RADAR MAST TO CABLE STANCHION IN FEET 
TO THE NEAREST FOOT. 
COL 10-16 ESTIMATED BACKGROUND COUNT RATE IN COUNTS/SEC FOR 
CHANNEL le 
COL 17-23 ESTIMATED BACKGROUND COUNT RATE IN COUNTS/SEC FOR 
CHANNEL 2e : 
COL 25-28 TIME OF INJECTION IN HOURS AND MINUTESe 
COL 29-31 THE NUMBER OF DAYS SINCE THE INJECTION 
COL 32-36 THE HALF LIFE OF THE RADIOISOTOPE IN DAYS. 
THE DECIMAL POJNT MUST BE PUNCHED INe 
COL 37-40 THE TIME WHEN THE SURVEY WAS STARTENe IN HOURS AND 
MINUTES 
COL 44-50 BEACONY NORTH COORD/ LAMBERT COORDINATES OF RADAR 
COL 52-58 BEACONI1 EAST COORD/ BEACONS TO THE NEAREST FOOT. 
COL 60-66 BEACON2 NORTH COORD/ BEACON IS ALWAYS UPCOASTeo 
COL 68-74 BEACOII2 EAST COORD/ 
COL 75-77 THE NUMRER OF LINES OF DATA TO BE SKIPPED AT THE 
BEGINNING OF A DATA SET IN ORDER TO AVOID 
READING IN SOME BAD DATAo 
FORMAT AND ENTRIES ON PLOT CONTROL CARD 
COL 1- 3. PLOT OPTION CONTROLe TO USE*e PUNCH THE NUMERAL 1 
IN THE COLIMN INDICATED. 
1 = PLOT TRACKLINE FOLLOWED BY SURVEY VESSELo 
2 = PLOT UNCORRECTED RADIATION VALUES 


3 PLOT RADIATION VALUES CORRECTED FOR BACKGROUND 
AND DECAY SINCE TIME ZEROe 
4 = UNUSED. LEAVE BLANKe 


COL 5-14 MAP SCALE EXPRESSED IN UNITS PER INCH 
COL 16-17 OPTION 1 / USE WHEN SPOTTING DATA FOR EACH PLOT 
COL 16-19 J 2 / OPTIONe USER CAN SPECIFY THAT EVERY NTH 
COL 20-21 J 4 / POINT BE PLOTTED. IF LEFT BLANKe THE 
PROGRAM ASSUMES EVERY POINT IS TO BE 
PLOTTED. 
COL 25-32 INTERVALS ON THE COORDINATE GRID AT WHICH TICK 
MARKS WITH THE LAMBERT COORDINATES WILL BE POSTEDe 
IF FIELD IS LEFT BLANKe PROGRAM WILL ASSUME THAT 
NO TICK MARKS ARE TO BE PLOTTED AND POSTED. 
COL 34-43 BEACON 1 / INJECTION SITEe DISTANCE IN METERS TO 
COL 45 54 BEACON 2 / THE NAMED BEACONSe IF FIELDS ARE LEFT 
BLANKe SITE IS NOT PLOTTED. 
COL 56-61 OAYeMONTH AND YEAR THE SANO WAS INJECTED 
FORMAT AND ENTRIES ON PLOT IDENTIFICATION CARD 
COL 1-78 THIS FIELD WILL BE PLOTTED ON THE LOWER MARGIN OF 
THE MAPe 
COL 80 PUNCH *T® HERE ON THE LAST DATA SET. 

COMMON /AA/ NOPT(4) oNPLT (3) e SCALE» GRID» DUMPNe DUMPE epNBARe EBAReLEGEN 
1D(13) e INDATE*LINE» BGCReSIGMAeNMAX eNMINe EMAX eo EMINe JOBENDeNENTRY oNCA 
2LLS»NCAM 

REAL NORKTH(2000) eNCORD(2000) »NBAReNMAX eNMIN 

DIMENSION NMBR(2000)* TIME(2000)» 0(202000)» EAST(2000)» ECORD(200 
10)» RAD(2000)» CCR(2000)» FATH(2000) 

EGUIVALENCE (RAD(1)0CCR(1))° (D101) eNCORD(1))0 (D0¢191001) PECORD(1 
1))o 

LOGICAL JOBEND 
DATA CORR1?CORR2/20902¢8/IFLAG/1H*/ 

WRITE (60610) 

NCAM=9 

NCALLS=0 

NCALLS=NCALLS+1 
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PROGRAM == RAPLOT III MODIFIED FOR IBM=7094 WITH SC=4060 CRT 


READ IN DATA CONTROL CARD 


READ (59620) CABLE» DEPTH? BOAT + BKGe ZHR»e ZMINe DAYSeHLIFE? SETIMEeRMINOG 
1SEC eBEACINe BEAC1E? BEAC2Ne BEAC2Ee ISKIP 


READ IN PLOT CONTROL PARAMETERS 


READ (50630) NOPT»SCALEe (NPLT(1) e I=103) »GRID*eDUMP1 » DUMP2 > INDATE 
DO 30 I=103 

IF (NPLT(I)) 20020030 

NPLT(T)=1 

CONTINUE 


READ IN PLOT LEGEND 


READ (50640) LEGEND»? JOBEND 
WRITE (60650) LEGENDe JOBEND 
WRITE (60660) BEACINe BEACiE* BEAC2Ne BEAC2E 


COMPUTE PROGRAM PARAMETERS FROM DATA CONTROL CARD ENTRIES 


SQDSTB=(BEAC2N=BEACIN) **2+ (BEAC2E=BEACIE) +22 
DISTB=SQRT(SQDSTB) 

WRITE (60670) SGOSTBeDISTB 
SINE=(BEAC2N-BEACIN) /DISTB 
COSINE=(BEAC2E=BEACIE) /DISTB 

WRITE (60680) SINE* COSINE 

WRITE (60690) ZHReZMINeSETIME*RMINe SECeDAYS 
ZHR=ZHR+ZMIN/60 © 

SETIME=SETIME+RMIN/606 
DELAY=SETIME+DAYS*246—-ZHR 

IF (HLIFE*GTe060) DECAY=ALOG(20)/(HLIFE*24 0) 
WRITE (60700) HLIFE*DECAY»DELAY 

WRITE (60710) CABLE»DEPTHe BOAT 
CABLE=BOAT+SQRT (CABLE **2=DEPTH*##2) 

WRITE (60720) CABLE 

WRITE (60730) (NOPT(I) eo T=103) 


COMPUTE COORDINATES OF THE INJECTION SITE FROM THE DISTANCES FROM 
THE BEACONS 


IF (DUMP1) 70°70°40 
DUMP1=(DUMP1+CORR1) #3228083 
OUMP2=(DUMP2+CORR2) #3228083 
DX1=(SQDSTB+DUMP 1*DUMP1-DUMP2*DUMP2) /(DISTB*2.) 
DY 1=DUMP1*DUMP1=-DX1*DX1 

IF (DY1) 50°50°60 

SITEN=-999999. 

WRITE (60740) 

GO TO 70 

OY1==-SQRT(DY1) 
SITEE=DX1*COSINE=DY1*SINE+BEACIE 
SITEN=0X1*SINE+DY1*COSINE+BEACIN 
WRITE (60750) SITEN*SITEE 


READ IN THE OATA FILE FROM ONE RIST SURVEYe ELIMINATING ANY LINES 
THAT CONTAIN A =30 WHICH IS AN ERROR FLAG 


N=1 
CALL RESET . 
READ (50760) NMBR(N) oe TIME(N) eD(L0N) 0D(20N) PRAD(N) oFATH(N) oe JFLAG 


CHECK FOR END OF FILE 


CALL CHECK (E) 
IF (E) 1006900100 


CHECK ERROR FLAG ON INPUT RECORD 


IF (IFLAGeNEoJFLAG) GO TO 80 
IF (NMBR(N) eLT20) GO TO 80 
IF (TIME(N) cL T0000) GO TO 80 
IF (D(1eN) oL T2020) GO TO 80 
IF (D(2eN) oLTe0.0) GO TO 80 
IF (RAD(N) oLTe000) GO TO 80 
IF (FATH(N) oLTe000) GO TO 80 
IF (NeGEe.000) GO TO 100 
N=N+1 

GO TO 80 
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PROGRAM == RAPLOT III MODIFIED FOR IBM~-7094 WITH SC-4060 CRT 


CONT INUE 

N=Nw1 

LINE=N-1 

IF (ISKIP-LE-0) GO TO 130 


SKIP LEADING CARD IMAGES THAT CONTAIN BAD DATAe 


NSTART=ISKIP#1 
DO 120 N=NSTARTeLINE 
NN=N=ISKIP 

NMBR (NN) =INMBR ON) 
TIME(NN) =TIME(N) 
DO 110 I=lee 
D(LTetN)=O0(TeN) 
RAD(NN) =RAD(N) 
FATH(NN) =FATH(N) 
LINE=LINE-ISKIP 
MSTOP=LINE=1 


CHECK DISTANCES TO RADAR BEACONS FOR ERRORSe IF DISTANCE/TIME 
FOR SUCCESSIVE BEACON RANGES INDICATE A SHIP SPEED eGTe 6 KNOTS 
(3009 METERS/SEC) PRANGE IS IN ERROR: 


bO 240 I=le2 

DO 170 M=lerSTOP 

IF (O(IoM)) 17001700140 
NSTART=M+1 

DO 160 N=NSTARToOLINE 

IF CABS(D(ITeN) -D( 10M) )=(TIME(N) ©TIME(M) ) #3208865) 17001700150 
O(LleN)=—1e 

CONTINUE 

CONTINUE 

00 230 M=1leMSTOP 

IF (O(IeM)) 18001807230 
NSTART=M 


CORRECT ERRONEOUS BEACON RANGES BY LINEAR INTERPOLATION (ON TIME) 
BETWEEN NON-ERRONEOUS RANGESe 


DO 200 N=NSTARTeLINE 

IF (O(IeN)) 200702000190 

NSTOP=N 

60. TO 210 

CONTINUE 

OTIME=TIME (NSTOP)-TIME(NSTART@=1) 
O0D1=D0(1eNSTOP)=D( I oNSTART=1) 
N=NSTART 

OC IeoN)=D( I oNSTART=1) +001 % (TIME (N) -TIME (NSTART=1) ) /OTIME 
N=N+1 

IF (N=NSTOP) 22002300230 
CONTINUE 

CONTINUE 


COMPUTE POSITION OF SHIP FROM DISTANCES FROM THE TWO BEACONS 


LAG=0 
DO 270 N=1leLINE 


MAKE CONSTANT CORRECTION FOR CUBIC AUTOTAPE INTERROGATOR 
AND CONVEKT TO FEET 


DFT1=(0(10t) +CORR1) #3228083 
OFT2=(D(2°N) +CORR2) #3.28083 
DX1=(SQUSTB+OFT1*OFT1-DFT2*DF T2) /(20*DISTB) 
DY1=DFT1*OFT1-DX1*DXxX1 


CHECK FOR IMAGINARY ROOTs 


IF (0Y1) 25002500260 
NORTH(N)=~le 
LAG=LAG+1 

GO TO 270 
OY1=-S@RT(DY1) 


ROTATE COORDINATES AND TRANSLATE INTO CALIFORNIA LAMBERT COORDINAT 
SYSTEM 


EAST(N)=DX1*COSINE-DY1*SINE+BEACLE 
NORTH(N) =DX1*SINE+DY1*COSINE+BEACIN 
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PROGRAM =—- RAPLOT III MODIFIED FOR IBM-7094 WITH SC=4060 CRT 


CONTINUE 
IF (LAGeEQe0) GO TO 330 
ASSIGN 330 TO KEY 


ELIMINATE DATA SETS FOR WHICH AN 
IMAGIN ARY FIX WAS OBTAINED 


LAG=0 

DO 320 N=leLINE 

IF (NORTH(N)) 29002900300 
LAG=LAG+t1 

GO TO 320 

NN=N=LAG 

NMBR (NN) =NMBR(N) 
TIME (NN) =TIME(N) 

DO 310 I=1+r2 
D(TeinNY=O0(IeN) 

RAD (NN) =RAD(N) 
EAST(NN) =EAST(N) 
NORTH(NN) =NORTHI(N) 
CONTINUE 
LINE=LINE*LAG 

GO TO KEYe (330+°400) 


CHECK NORTI4 AND EAST COORDINATES FOR 
EXTREME VALUES BY CHAUVENET'S CRITERION 


NBAR==-1. 
CALL STDEV (NORTHeLINE*NBAR® SONRTH) 
EBAR=-1. 


CALL STDEV (EASTeLINE*EBARe SDEAST) 
ALPHA=1.0-1-/FLOAT (2*LINE) 

CALL TINORM (CHVR»e ALPHA? $335) 

GO TO 340 

CHVR=5.6 

WRITE (60770) ALPHA 
GATE1=EBAR=CHVR*SDEAST 
GATE2=EdAR*+CHVR*SDEAST 
GATEN1=NBAR=CHVR*SONRTH 
GATEW2=NBAR+CHVR*SDNRTH 

LAG=0 

DO 390 W=leLINE 

IF (EAST(N)-GATEL) 380°350°350 
IF (EAST(N)=GATE2) 360°360+°380 
IF (NORTH(N)=GATEN1) 38003700370 
IF (NORTH(N)=GATEN2) 39003900380 
NORTH(N)=-1e 

LAG=LAG+1 

CONT INUE 

ASSIGN 400 TO KEY 


ELIMINATE ANY DATA SETS THAT HAVE AN EXTREME 
VALUES OF THE NORTH OR EAST COORDINATES 


IF (LAGeGTe0) GO TO 280 
CALL SUBROUTINE FOR PLOTTING THE TRACK OF THE SURVEY VESSELe 


NENTRY=1 
IF (NOPT(1)EQ@.1) CALL TRACK (NORTHEAST oNMBR) 


APPLY A CORRECTION TO ALLOW FOR THE DISTANCE THE DETECTOR IS TOWED 
ASTERN OF THE SURVEY SHIP. 


DNO=NORTH(1)=(NORTH(2) =NORTH(1) ) 
DEO=EAST(1)-(EAST(2)-EAST(1)) 

DENOM=SQRT ( (NORTH(1)=DNO) «#24 (EAST(1)-DEO) *#2) 
NCORD(1)=NORTH(1)-CABLE* (NORTH(1)—-DNO) /DENOM 

ECORD(1) =EAST (1) -CABLE* (EAST (1)-DEO) /DENOM 

DO 430 N=2eLINE 

DENOM=SORT ( (NORTH(N) =NCORD(N=1) ) ##2+4 (EAST(N) -ECORD(N=1) ) #2) 


THE CORRECTION FOR THE DISTANCE BETWEEN VESSEL AND THE DETECTOR 
IS EQUAL TO "CABLE® UNLESS THE VESSEL IS LESS THAN °CABLE® FEET 
AWAY FROM THE LAST COMPUTED POSITION OF THE DETECTOR VEHICLE. IN 
THIS EVENTe THE NEW COMPUTED DETECTOR POSITION IS THE SAME AS 

THE LAST DETECTOR POSITION 
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PROGRAM -- RAPLOT III MODIFIED FOR IBM-7094 WITH SC~4060 CRT 


IF (DENOM=CABLE) 420¢420°410 

NCORD (N) =NORTH(N) =CABLE* (NORTH(N) -NCORD(N=1) ) /DENOM 
ECORD(N) =EAST(N) -CABLE* (EAST (N) -ECORD(N=1) ) /DENOM 
GO TO 430 B 
NCORD(N)=NCORD(N=1) 

ECORD(N)=ECORD(N=-1) 

CONT INUE 


CONVERT RADIATION READINGS TO COUNTS PER SECOND 


CCR(1)=RAD(1) /SEC 

DO 450 N=2eLINE ; 

IF (NMBRON) eNEeNMBR(N=1)41) GO TO 440 
CCR(N) =RAD(N) / (TIME (N) <TIME (N=1)) 

GO TO 450 

CCR(N) =RAD(N) /SEC 

CONTINUE 


COMPUTE THE MEAN AND STANDARD DEVIATION OF THE BACKGROUND COUNT 
RATE FROM THE RADIATION DATA THAT LIES WITHIN THE LIMITS OF THE 
ESTIMATED BACKGROUND COUNT RATE SET BY CHAUVENET’S CRITERION. 


IF (CCR(1)eGTe2e*BKG) CCR(1)=BKG 
ALPHA=1.—-10/FLOAT(2*LINE) 
CALL TINORM (CHVRe ALPHA?$455) 
GO TO 460 

CHVR=5. 

WRITE (60770) ALPHA 

BGCR=0.0 

SIGMA=0.0 

NCOUNT=0 

GATE1=BKG=CHVR*SGRT (BKG/SEC) 
GATE2=BKG+CHVR*SQRT (BKG/SEC? 
DO 490 N=1eLINE 

IF (CCR(N)=GATE1) 49004900470 
IF (CCR(N)=GATE2) 48004909490 
NCOUNT=INCOUNT+1 
BGCR=BGCRtCCRIN) 

CONTINUE 

BGCR=BGCR/FLOAT (NCOUNT) 

DO 520 N=1leLINE 

IF (CCR(N)=GATE1) 520°520/500 
IF (CCR(N)-GATE2) 5100520520 
SIGMA=SIGMA+(CCR(N) =BGCR) ¥#2 
CONTINUE 

SIGMA=SGRT (SIGMA/FLOAT (NCOUNT) ) 
WRITE (60780) BKGeBGCReSIGMA 


CALL THE SUBROUTINE FOR PLOTTING UNCORRECTED RADIATION VALUES 


NENTRY=2 
IF (NOPT(2)eE@e1) CALL RADPLT (NCORDeECORDeCCR) 
IF (NOPT(3) eNEe1) GO TO 570 


CORRECT RADIATION VALUES FOR BACKGROUND COUNT RATE AND TIME=DECAY 


SUM=0.0 

NBAR=0.0 

EBAR=0-0 

SDNRTH=0 0 

SDEAST=0.0 

DO 540 N=leLINE 

CCR(iV) =CCR(N) -BGCR 
IF(CCR(N)=30e*SIGMA) 54005400530 
CCR(N)=(CCRIN) —30*SIGMA) EXP (DECAY® (DELAY+TIME(N)/36000) ) 
1 + 36¥*SIGMA 

SUM=SUM+CCRI(N) 
NBAR=NBAR+ (NCORD(N) =NCORD(1))#CCR(N) 
EBAR=EBAR+ (ECORD(N) -ECORD(1)) #CCRIN) 
CONTINUE 

WRITE (60790) SUM 


COMPUTE WEIGHTED MEAN AND STDe DEVe OF ACTIVITY LOCA) suN 
NBAR=NCORD (1) +N6BAR/SUM 
EBAR=ECORD(1)+EBAR/SUM 


DO 560 N=1eLINE : 
IF (CCR(N)) 560°560°550 
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PROGRAM == RAPLOT III MODIFIED FOR IBM=-7094 WITH SC-4060 CRT 


SDNRTH=SDNRTH+ (NCORD(N) -NBAR) * (NCORD(N) -NBAR) #CCR(N) 
SDEAST=SDEAST+ (ECORD(N) -EBAR) * (ECORD(N) -EBAR) #CCR(N) 
CONTINUE 

SONRTH=SQRT (SONRTH/SUM) 

SDEAST=SQRT (SDEAST/SUM) 

WRITE (60800) NBAR*EBARs SONRTH?e SDEAST 


COMPUTE AND PRINT 95 PCe CONFIDENCE LIMITS OF MEAN RADIATION 
LOCATION. 


RTSUM=SQRT (SUM/BGCR) 

CF IDN=1 e96*SONRTH/RTSUM 

CF IDE=1.96*SDEAST/RTSUM 

WRITE (60810) CFIDN*eCFIDE 

NENTRY=3 

IF (NOPT(3)eEQ@e1) CALL RADPLT (NCORDeECORDe CCR) 
WRITE (60820) NMAX*®EMAX 

WRITE (60830) NMINeEMIN 


WRITE OUT THE NUMBERe COORDINATES AND ACTIVITY OF EACH DATA POINT 


KOUNT=50 

DO 600 N=leLINE 

IF (KOUNT=50) 59005800580 

WRITE (60840) LEGEND 

KOUNT=0 2 
WRITE (69850) NMBR(N) oe TIME(N) eNORTH(N) oEAST(N) eNCORD(N) PECORDIN) oC 
ACR(N) ePFATH(N) 

KOUNT=KOUNT+1 

WRITE (60860) 

IF (.NOTeJOBEND) GO TO 10 

STOP 


FORMAT (1H1) 

FORMAT (3F3c00F7e00 8X0 2F2.09F3c0PFS5ce20 3F20004(1X0F 700) 013) 

FORMAT (4I1eF10e0e1XeST201X#F10eOr2(1X0F 1000) 0 1X0 AG) 

FORMAT (13A6eL2) 

FORMAT (10X013A6010X0L2) 

FORMAT (/5Xe8HBEACON 1eF10¢001HNeF10¢0»LHEeSXeBHREACON 20F10.0°1HN 
LeF1l0eQvelHE//) ~~ ‘ 

FORMAT (5Xe21HSQUARE DIST BETWEEN =0E16.895Xe1BHDISTANCE BETWEEN = 
1°F10-0) 

FORMAT (SX*6HSINE =cE£16¢8/5XeBHCOSINE =0£1608) 

FORMAT (5Xe16HIWJECTION TIME =02F3.0°10Xe12HCLOCK SET AT02F3.005Xe 
121HDIGITIZING INTERVAL =0F3e0* THSECONDS/5Xe30HDAYS ELAPSED SINCE I 
2NJECTION =0F3-0) 

FORMAT (SX*e22HHALF=LIFE OF ISOTOPE =eF7e2°4HDAYS*?5Xe 14YHDECAY FACTO 
1R =0&16¢8°S5Xe19HTIME=DELAY FACTOR =0F7.205HHOURS) 

FORMAT (10Xe14HCABLE LENGTH =0F4e0010Xe18HMEAN WATER DEPTH =0F4.00 
110X°13HBOAT LENGTH =°F4e0) 

FORMAT (10X¢38HDISTANCE FROM RADAR MAST TO DETECTOR =eF6.10SHFEETe 
1) 

FORMAT (//40Xe15HPLOTS GENERATED/10Xe 9SHTRACKLINE oe 150 10X » 20HBACKGRO 
LUND RADIATION? 15010Xe19HCORRECTED RADIATION®e I5) 

FORMAT (SXeS1HBEACON RANGES FOR DUMP SITE COMPUTE IMAGINARY ROOT.) 

FORMAT (10X»37HLAMBERT COORDINATES OF INJECTION SITEeF10.001HNoF10O 
1-0°1HE) 

FORMAT (I603F7cel0F 7000 TXOFT COP AL) 

FORMAT (37H THERE WAS AN OVERFLOW WHEN ALPHA WASeF6e3e5Xe25HCHVR W 
1AS SET EQUAL TO 5-0) 

FORMAT (//20X*S3HSUMMARY STATISTICS OF BACKGROUND RADIATION COUNT 
URATE/10Xe21HEST* BKGe COUNTS/SECeeF10¢0/10Xe21HMEAN BKGe COUNTS/SE 
2CorF10e0/10X021HSTD» DEVe COUNTS/SEC. F100) 

FORMAT (//10Xe33HSUM OF CORRECTED RADIATION COUNTSeE16-8) 

FORMAT (//20Xe41HSUMMARY STATISTICS OF RADIATION LOCATION. /24X011H 
UNORTH COORD? LOX» LOHEAST COORD/16Xe4HMEAN® SX0F 10-00 10X0F1000/11X09H 
25TDe DEVe re SXeF10c¢Or01L0XeF 10.0) 

FORMAT (10X° LOHCONFIDENCE/7X013HLIMIT OF MEAN®5X0F10-0010X0F 1000) 

FORMAT (//7Xe13HMAXIMUM COORD¢5X0F 10-00 1HNe 9X0F 10000 LHE) 

FORMAT (//7Xe0L3HMINIMUM COORD¢5XeF10¢00 LHNe9X0F10¢09 LHE/S/) 

FORMAT (1H1°9X013A6//2XeS5OHLINE TIME DISTANCE TO BEACON BOAT 
ACOORDTNATES*?4X»SOHBALL COORDINATES UNCORRECTED CORRECTED DEP 
2TH/9K» SHSEC 8X0 LHL 0 9X0 1LH20 6X 0 SHNORTH? 6X eo 4HEAST 0 5X o SHNORTH9 6X9 4HEAS 
3T°4Xe20HRADIATION COUNTS/SEC 0 6X0 4HFEET) 

FORMAT (1X9 150F6c0020X04F1000013X0F11-00F10-0) 

FORMAT (1H1) 

END 
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467 
468- 


10 

20 

30 

40 

50 

60 

70 

80 

90 

100 
110 
120 
130 
140 
150 
160 
170 
180 
190 
200 
210 
220 
230 
240 
250 
260 
270 
280 
290 
300 
310 
320 
330 
335 
340 
350 
360 
370 
380 
390 
400 
410 
420 
430 
440 
450 
455 
460 
470 
480 
490 
500 
510 
520 
530 
540 
550 
560 
570 
580 
590 
600 
610 
620 
630 
640 
650 
660 
670 
680 
690 
700 
710 
720 
730 
740 
750 
76U 
770 


PROGRAM == RAPLOT III MODIFIED FOR IBM=7094 WITH SC-4060. CRT 


10# 
14 
13 
36 
41 
41 
36 
50* 
53 
53 
73 
69 
67 
61 
84 
83 
80 
89 
92 
91 
94 
99% 
88 
79 
110 
110 
105 
120% 
122 
122 
128 
121i 
118 
143% 
142 
151 
152 
153 
151 
150 
135 
169 
169 
167 
178 
177 
187% 
186 
195 
196 
194 
202 
203 
201 
218 
216 
228 
227 
210 
245 
245 
244 
7WR 
11RD 
12RD 
17RO 
18WR 
19WR 
22WR 
25wR 
26wWwR 
31wR 
32WR 
34WR 
35WR 
43WR 
4BWR 
51RD 
144WR 


251i 
15* 
14 
37% 
4axe 
45% 
4y 
54 
54x 
61 
74x 
76% 
78x 
82% 
65% 
86* 
81 
90% 
93% 
92 
96% 
101 
89 
103% 
11i1x* 
114% 
113 
159 
123% 
125% 
129% 
124 
119 


145% 
152% 
153% 
154% 
152 
154 
158 
170* 
173% 
172 
181% 
180 


189% 
196* 
197% 
195 
203% 
204% 
202 
219% 
218 
229% 
228 
239% 
246% 
248% 
249" 
253* 
254% 
255* 
256* 
257% 
258* 
259* 
260* 
261* 
262% 
263% 
264% 
265* 
266% 
267* 
268% 
1Aa8WR 


16* 


49x 
s5 


64" 


64 


95% 


117% 


133% 
135 


153 

157% 
160* 
175% 


182% 


196 


203 
223% 


231% 


269% 


56 57 58 59 60 


87% 


102% 


136% 


154 155% 


199% 


205% 
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PROGRAM == RAPLOT III MODIFIED FOR IBM=-7094 WITH SC=8060 CRT 


780 = 207WR 270% 

790 = 224WwR 271* 

800 = 234WR 272s 

B10 - 238WR 273% 

820 = 241WR 274 

830 = 242wR 275* 

840 - 246WR 276 

850 = 24aWR 277% 

860 = 250WR 278s 

AA - 1c0 

ABS - 84 

ALOG - 30 

ALPHA - 140= 141AG  144WwR =184= 185AG  188WR 

BEACLE = 11RD 19WR 20 au 46 115 

BEACIN = 11RD 19WR 20 23 47 116 

BEAC2E = 11RD 19WR 20 24 

BEAC2N = 11RD 19WR 20 23 

BGCR - 1cO. =: 189= 198= 2002 204 207WR = 217 235 

BKG - 11RD 183 192 193 207WR 

BOAT - 11RD 32WR 33 

CABLE = 11RD 32wR 33= 34WR = 165 166 169 170 171 

CCR - 301 4EQ 176= 179= 181= 183 195 196 198 202 
203 204 209AG =. 217= 218 219= 220 221 222 228 
229 230 240AG 246WR 

CFIDE <- 237= 238wR 

CFION - 236= 238WR 

CHECK = 52 

CHVR - 141AG 143= 145 146 147 148 165aG =: 187= 192 193 

CORRL = 6DA 37 106 

CORR2 - 6DA 38 107 

COSINE - 24u= 25wR 46 47 115 116 

D - 301 4EQ 51RD 57 58 742 81 84 85= 89 
92 97 99= 106 107 129= 

DAYS - 11RD 26W0R 29 

DD1 - 97= 99 

DEO - 163= 164 166 

DECAY - 30= 3lwR 219 

DELAY = 29= 3iwR 219 

DENOM = 1642 165 166 168= 169 170 i71 

DEPTH = 11RD 32wR 33 

OFT1 -  106= 108 109 

DFT2 -  107= 108 

DISTR. - 21= 22WR 23 a4 39 108 

DNO - 162= 164 165 

DTIME - 96= 99 

OUMPL = - 12RD 36 37= 39 40 

DUuMP2 = 12RD 3a= 39 

DUMPE - 1C0 

DUMPN = 1co 

Dx1 - 39= 40 46 47 108= 109 115 116 

pY1 - 4o= 41 45= 46 47 109= 110 114= 115 116 

E - 52AG 53 

EAST - 301 115= 131= 139AG 6:15 152 161AG 163 164 166 
168 171 248WR 

EBAR - ico. )—-:138= 139aG 145 146 213= 222= 226= 230 234WR 

EcoRD - 301 4EQ 166= 168 17i= 174= 2090G «= 222 226 230 

240AG = 248wR 

EMAX - ico. =. 241 wR 

EMIN - ico 242WR 

EXP - 219 

FATH - 301 51RD 60 76= 24B8WR 

FLOAT - 140 184 200 206 

GATEL - 145= 151 192= 195 202 

GATE2 = 146= 152 193= 196 203 

GATEN] = 147= 153 

GATEN2 = 148= 154 

GRID - 1co 12RD 

HLIFE = 11RO 30 31WR 

I - 12R0 13 14 15 35WR 73 74 79 81 a4 

85 89 92 97 99 128 129 

IFLAG - 6DA 54 

INDATE = 1co 12RD 

ISKIP = 11RD 67 68 70 77 

JFLAG = 51RD 54 

JOBEND = 1co 5LG 17RD 18WR 251 

KEY - 119= 135 158= 

KOUNT = 243= 245 247= 249= 

LAG = 104= 112= 118 120= 123= 125 134 149= 156= 159 
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PROGRAM == RAPLOT III MODIFIED FOR IBM-7094 WITH SC=4060 CRT 


LEGEND = 1co 17RD 18wR 246wWR 

LINE = ico 66= 69 77= 78 83 U1 105 121 134= 
137AG 139AG 150 167 177 194 201 216 227 244 

M = 80 61 82 84 88 69 90 

MSTOP 2 78= 80 88 

N - 49= 51RD 55 56 57 58 59 60 61 62= 
65= 66 69 70 71 72 74 75 76 83 
84 85 91 92 93 98= 99 100= 101 105 
106 107 111 115 116 121 122 125 126 127 
129 130 131 132 150 151 152 153 154 155 
167 168 170 171 173 174 177 178 179 181 
194 195 196 198 201 202 203 204 216 217 
218 219 220 221 222 227 228 229 230 244 
248WR 

NBAR = 1co 2RL 136= 137AG 147 148 212= 221= 225= 229 
234WR 

NCALLS = 1co 9= 10= 

NCAM aa 1co = 

NCORD 2 2RL 4EQ 165= 168 170= 173= 209AG 221 225 229 
240AG 248wR 

NCOUNT = 1g91= 197= 200 206 

NENTRY = 1co 160= 208= 239= 

NMAX = 1co 2Rut 241wWR 

NMBR a 3DI 51RD 55 71= 126= 161AG 178 248WR 

NMIN O 1co 2RL 242WR 

NN - 70= 71 72 74 75 76 125s 126 127 129 
130 131 132 

NOPT 2 1co 12RD 35WR 161 209 210 240 

NORTH = 2RL 111= 116= 122 132= 137AG 153 154 155= 161AG 
162 164 165 168 170 248WR 

NPLT © 1co 12RD 14 15= 

NSTART = 68= 69 62= 83 90= 91 96 97 98 99 

NSTOP Oo 93= 96 97 101 

RAD oO 301 4EQ 51RD 59 75= 130= 176 179 181 

RADPLT = 209 240 

RESET = 50 

RMIN = 11RD 26WR 28 

RTSUM - 235= 236 237 

SCALE oO ico 12RD 

SDEAST = 139AG 145 146 215= 230= 233= 234wWR 237 

SDNRTH .= 137AG 147 198 2145 2295 2325 234WR 236 

SEC 2 11RD 26WR 176 161 192 193 

SETIME = 11RD 26wWR 28= 29 

SIGMA - ico 190= 2045 206= 207Wk 218 219 

SINE = 23= 25wWR 46 47 115 116 

SITEE a 46= 4BWR 

SITEN Oo 42s 47= 48WR 

S@DSTB = 20= 2l 22WR 32 108 

SORT = 21 33 4s 114 164 168 192 193 206 232 
233 235 

STDEV = 137 139 

STOP = 252 

SUM = 211= 220= 224WR 225 226 232 233 235 

TIME O 3DI 51RD 56 72= 64 96 99 127= 179 219 
246WR 

TINORM = 141 185 

TRACK a 161 

ZHR = 11RD 26WR 27= 29 

ZMIN o 11RD 26WR 27 


mefotetotoepotpatetoatapetotoetopotototoepototopofpotutopototetofpopotopatototatetofotoatoteopototoy 
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fe 
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50 


60 


70 


aan @a2nan 


aan 


100 
110 


ana00 


SUBROUTINE TRACK (NORTHeEASTeNMBR) 


SUBROUTINE TRACK (NORTH*EASTeNMBR) 


THIS SUBROUTINE GENERATES PLOT INSTRUCTIONS FOR THE STROMBERG- 


CARLSON 4020 CATHODE RAY TUBEe 


TRK 
TRK 
TRK 
TRK 
TRK 


COMMON /AA/ NOPT(4) eNPLT (3) oe SCALE’ GRIDe DUMPNe DUMPE e NRAR eEBAR e LEGENTRK 
1D (13) e INDATE*NPTSeBGCRe SIGMA’ NMAX e NMINe EMAX e EMINe JOBENDeNENTRY eNCATRK 


2LLSeNCAM 

DIMENSION EAST(2000)» NMBR(2000) 

REAL NORTH(2000) eNBAReNMAX eNMINeNMAXV 
LOGICAL JOBEND 

IF (NCALLSeGTe1l) GO TO 20 


INITIALIZE CAMERAe AND PRINT ID FRAME 


PHI=ALOG(2.) 

CALL FRAMEV 

CALL SETMIV (0014003200) 
LAST=3 

IF (NOPT(LAST) eNE2O) GO TO 40 
LAST=LAST~1 

GO TO 30 


COMPUTE MAXIMUM AND MINIMUM VALUES OF NORTH AND EAST ARRAYS 


NMAX=NORTH(1) 

NMIN=NORTH(1) 

EMAX=EAST(1) 

EMIN=EAST (1) 

DO 50 N=2eNPTS 
NMAX=AMAX1 (NMAX e NORTH(N) ) 
NMIN=AMIN1 (NMINe NORTH(N) ) 
EMAX=AMAX1 (EMAXe EAST (N) ) 
EMIN=AMIN1(EMINeEAST(N) ) 
XRATIO=(EMAX=EMIN) /8830 
YRATIO=(NMAX=NMIN)/991 0 

IF (XRATIO=YRATIO) 60980070 
EMAXV=EMIN+YRATIO*883. 
NMAXV=NMAX 

GO TO 6BU 
NMAXVENMIN4#XRATIO#991 © 
EMAXV=EMAX 

GO TO 80 

ENTRY RADPLT (NORTH? EAST RAD) 
DIMENSION RAD(2000) 

IF (NOPT(1) eEQeOeANDeNENTRYcEQe2) GO TO 10 
IF (NOPT(1) cE@e0eANDoNOPT (2) eEQ@eOeANDeNENTRY-EQ.3) GO TO 10 


PRINT LEGEND AND INJECTION DATE AT BASE OF MAPo 


CALL PRINTV (78eLEGEND? 24224) 

CALL PRINTV (-14°1GHINJECTION DATE? 72008) 

CALL PRINTV (6* INDATE*? 84808) 

NN=NPLT (NENTRY) 

GO TO (90¢110°140)» NENTRY 

CALL PRINTV (-25°25HPLOT OF SURVEY TRACK LINE®2408) 


CALL GRIDIV (2, EMIN’e EMAXV eNMINe NMAXVe SCALE 0 SCALE 05 0S0SeSe-Se-4) 


PLOT TRACK LINE FOLLOWED BY SURVEY VESSEL 
DO 100 N=2eNPTS 


TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 


CALL LINEV (NXV(EAST(N=1)) eNYV(NORTH(N=1) ) oNXV(EAST(N) ) oNYV(NORTH(TRK 


1n))) 

IF (MOD(NeNN) NEO) GO TO 1U0 
CALL POINTV (EAST(N) eNORTH(N) 2-3) 
D=FLOAT (NMBR(N) ) 

IX=NXV (EASTON) +5 


PLOT FIX NUMBER OF EVERY NNTH FIX 
CALL LABLV (De IXeNYV(NORTH(N) ) 040104) 


CONTINUE 
GO TO 220 


TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 


CALL PRINTV (-62°62HPLOT OF UNCORRECTED RADIATION IN STD. DEVSe FRTRK 


10M MEAN BKGe CRer24e8) 


CALL GRIDLV (2eEMINe EMAXV oNMINeNMAXVe SCALE 0 SCALE 050S5eSeSe—-Se—4) 


TRK 
TRK 
TRK 


PLOT UNCORRECTED RADIATION VALUES IN STANDARD DEVIATIONS FROM MEANTRK 


BACKGROUND COUNT RATE 
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TRK 
TRK 


ODNDNF UNE 


120 
c 


200 


ann 


SUBROUTINE TRACK (NORTHe EAST» NMBR) 


DO 120 N=leNPTS 

IF (MOD(NeNN) eNE*O) GO TO 120 
IVAL=IFIX((RAD(N) -BGCR) /SIGMA+60) 
IF (IVALeLTeO) IVAL=0 

IF (IVALeGTe1l1) IVAL=11 

IVAL=-IVAL 

CALL POINTV (EAST(N) pNORTH(N) e IVAL) 
CONTINUE 


PRINT LEGEND FOR SYMBOLS REPRESENTING UNCORRECTED VALUES 


CALL PRINTV (-6e6HLEGEND? 9120800) 

CALL POINTV (887°768°0° ANY) 

CALL PRINTV (-8e8HCCR eLE. 09030768) 

D=BGCR=5-*SIGMA 

CALL LABLV (009750 768060106) 

1Y=752 

DO 130 I=1ell 

N=-I 

CALL POINTV (8870IYeNe ANY) ; 

CALL PRINTV (-8e8HCCR eGT.e9030TY) 

CALL LABLV (De9750IYe60106) 

D=D+SIGMA 

IV=1Y-16 

GO TO 220 

CALL PRINTV (=38e38HPLOT OF CORRECTED RADIATION COUNT RATE*24e8) 
CALL GRIDIV (2eEMINe EMAXV eNMINe NMAXV 0 SCALE 0 SCALE 0S 05eSe5e-S0—-4) 


PLOT CORRECTED RADIATION COUNT RATE AS BACKGROUND IF oLTe 3 9TD. 
DEVSe FROM BACKGROUNDe IF COUNT RATE IS MORE THAN 3 STDe DEVS. 
BELOW BACKGROUNDe THE VALUE IS NOT PLOTTED AT ALLeo 


DO 190 N=leNPTS 

IF (MOD(NeNN) eNE°O) GO TO 190 

IF (RAD(N)+3¢#SIGMA) 19001600150 
IF (RAD(N)=3e¢*SIGMA) 16001600170 
IVAL=0 

GO TO 180 


IF THE COUNT RATE IS eG1e BACKGROUNDe THE VALUE IS PLOTTED ON A 
POWER OF 2 * 25 SCALE 


IVAL=IFIX(ALOG(RAD(N) 7100.) /PHI+30) 
IF (IVALeLTol) IVAL=1 

IF (IVALeGTo12) IVAL=12 

IVAL=-IVAL 

CALL POINTV (EAST(N) eNORTHIN) o IVAL) 
CONTINUE 


PRINT LEGEND FOR SYMBOLS REPRESENTING CORRECTED VALUES 


CALL PRINTV (=6*°6HLEGEND? 9120800) 
CALL POINTV (887° 768°0eANY) 

CALL PRINTV (°8°8HCCR oLE. 09030768) 
D=3-*SIGMA 

CALL LABLV (009750 768060196) 

D=256 

1Y=752 

DO 2v0 I=1012 

N=-I 

CALL POINTV (887°1YeNe ANY) 

CALL PRINTV (=8¢8HCCR eGEer9030l1Y) 
CALL LABLV (De975e01T¥e60106) 

D=D*26 

IY=1Y-16 

PLOT MEAN RADIATION LOCATION 
IX=NXV(EBAR) 

TY=NYV (NBAR) A 

CALL POINTV (IXeITYoOe ANY) 

CALL PRINTV (=-8e8HO RADBAReIXe LY) 
PLOT POSITION OF INJECTION SITE 

IF (DUMPNeGTeNMAXeOReDUMPNeLTeNMIN) GO TO 210 


IX=NXV (OUMPE) 
IY=NYV (DUMPN) 
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TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 
TRK 


111 
112 
113 
114 
115 
116 
117 


SYMBOL 


10 
20 
30 
40 
50 
60 
70 
80 
90 
100 
110 
120 
130 
140 
150 
160 
170 
180 
190 
200 
210 
220 
AA 
ALOG 
AMAX1 
AMIN 
ANY 
BGCR 
0) 


DUMPE 
DUMPN 
EAST 


EBAR 
EMAX 
EMAXV 
EMIN 
FLOAT 
FRAMEV 
GRID 
GRIO1V 
I 

IFIX 
INDATE 
IVAL 


Ix 
IY 


JOBEND 
LABLV 
LAST 
LEGEND 
LINEV 
MOD 


210 CONTINUE 
220° CALL FRAMEV (2) 
RETURN 
ENO 
- 6* 34 35 
- 6 10* 
- 11* 13 
~ 11 14% 
- 18 22% 
2 25 26% 
= 25 29% 
i 25 28 31 
- 40 41s 
- 43 45 50* 
- 40 52% 
- 54 55 61s 
= 68 74% 
- 40 76* 
= 80 61* 
= 60 61 82% 
- 61 84x 
= 83 68* 
- 78 79 80 
- 97 103% 
= 108 113 
= Si 75 114% 
= 2co 
- i 84 
- 19 21 
- 20 22 
- TOAG 99AG 
2 2co 56 65 
- 47= 49AG 65= 
102= 
a 2co 109 
co 2co 108 110 
- 1AG 301 16 
60AG 88AG 
= 2co 104 
- 2co 16= 2i= 
a 26= 30= G2AG 
= 2co 17= 225 
- 47 
- 8 114 
2 2co 
- 42 53 77 
= 68 69 97 
- 56 84 
2 2co 38AG 
= 56= 57 +58 
8BAG 
© 48= 49AG 1045 
> 67= 7O0AG 71AG 
105= 106AG 107AG 
2 2co SLG 115 
— 49 66 72 
— 10= il l2= 
= 2co 36AG 
2 4y 
2 4s 55 79 
a 18 19 20 
48 49AG 54 
80 61 84 
= 2co SRL 105 
— 2co 6 
— 2co 
= 2co 34 35 
= 2co GRE 14= 
— 4RLE 27= 29= 
= 1AG 301 47 


CALL POINTV (IXeIYoOe ANY) 
CALL PRINTV (-606HO OUMPrIXe IY) 


IF (CeNOTeJOBENDeOReNENTRY eNEoLAST) RETURN 


SUBROUTINE TRACK (NORTHeEASToNMBR) 


SUBROUTINE TRACK (NORTH®EASToNMBR) 


36% 


89% 


66AG 


17 


23 


S3AG 


23 


98 


596 


106AG 


72AG 


110= 


94 


115 


al 
55 
88AG 


39 
19= 
42AG 


REFERENCES 
72A6-73= 
21. 22 
30 
7706 
26 Q2AG 
60AG = 82= 
10706 =: 109= 
74= 96= 

T11AG 112AG 

101 
22 43 
56 60AG 
98= 9946 
40 115 
24 27 
53A6 776 


59 


93= 


32 


S3AG 


64s 


111AG 
99AG 


4UAG 
69= 


108 


94AG 


GUAG 


T7AG 


-85 


112AG 
100AG 


as 
70AG 


95= 


46AG 


86 


101AG 


G6AG 
78 


158 
159 
160. 
161 
162 
163 
164=- 


101AG 


48 


a7s 


103= 


47 
79 


SUBROUTINE TRACK (NORTHeEASTeNMBR) 


NMIN 2 2co SRL 15= 20= 24 29 &2AG 53AG T7AG 108 

NN 2 39= 45S 55 79 

NOPT 2 2co il 34 35 

NORTH = 1AG GRE 14 15 \9 20 32 S4AG 46AG6 B9AG 
60AG 88AG 

NPLT — 2co 39 

NPTS = eco 18 43 54 78 

NXV QO 44UAG 48 104 109 

NYV 2 44UAG 49AG 105 110 

PHI = 7= 84 

POINTV = 46 60 63 70 88 9o1 99 106 111 

PRINTV = 36 37 38 41 52 62 64 71 76 90 
92 100 107 112 

RAD > 32 3301 56 80 61 84 

RADPLT = 32 

RETURN = 115 116 

SCALE 2 2co 42AG 5S3AaG T7AG 

SETMIV = 9 

SIGMA — 2co 56 65 - 73 80 61 93 

TRACK = i 

XRATIO = 23= 25 29 

YRATIO = 245 25 26 


efofaototopofototofpopsopofopoypotosofopotofopopotopofosotafopofpofpopopopotpopopgopotop¢ofopopogoy 


INDEX 
SUBROUTINE TINORM(Z2VALe ALPHA: *) 
1 SUBROUTINE TINORM(ZVAL?e ALPHA? #) TINORM 
2 DIMENSION A(3)0 B(3) TINORM 
3 DATACACL) oe 1=103) 700103289 080285302.515517/0 (B(1) of=103)/.00103089 TINORM 
1.189269010432788/ TINORM 
(Pete rn een nn nnn n non ee en eres eset esate sere sess eseesteresesetasesesete2=<7 INORM 
c APPROXIMATION TO INVERSE NORMAL DISTRIBUTION TINORM 
Ceeteesecerese===— eee meme en eee eee ern eee meen eee oe eee eee ee ee2eee2ece2e222T I NORM 
4 IF (eNOTo (ALPHA cGT 000 cANDcALPHAcLTele)) GO TO i TINORM 
5 X=ALPHA TINORM 
6 IF (XeGT.oe5) X=1Le-X TINORM 
7 X=SORT (2 ¢ *ALOG(X)) TINORM 
8 ZVAL =X=(A(3)4X#(A(2)4X2A(1)) 7 (10 4X8 (B13) 4xXe(B(2)4X* BC1)))) TINORM 
9 CALL OVERFL(1) TINORM 
10 IF(1eEQ@.1) RETURN 3 TINORM 
il IF(ALPHAcLTeo5) ZVAL = - Z2VAL TINORM 
12 RETURN TINORM 
13 1 RETURN 3 TINORM 
14 END TINORM 
INDEX 
SUBROUTINE TINORM(ZVAL¢ ALPHA: *) 
SYMBOL SPeseseeueBeesesessesea REFERENCES sSs2S2e2222222235228 
1 = 4 13% 
A o 201 3DA 8 
ALOG = U 
ALPHA ) 1AG 4 5 il 
B = 201 3DA 8 
1m - 3DA 9AG 10 
OVERFL - 9 é 
RETURN = 10 12 13 
S@RT o 7 
TINORM = 1 
x = = 6 7= 8 
ZVAL = 1AG 6= 11= 


epedatatetetatetepapotototoegetototetotetataototototetotetatqopaotatatatatotatototetotototateat 
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consists of ranges to two shore-based radar beacons, are first edited for spurious 
data, and then converted to rectangular coordinates (in this case the California 
Lambert Coordinate System). Radiation data are converted to count rate as counts 
per second. Background count rate is computed and subtracted from the observed 
count rate, and any radiation counts that are significantly above the background 
count rate are corrected for time of decay since the isotope was injected. 


Output from the program is in three forms - printed output, graphical output, 
and magnetic tape record. The processed data are transferred to magnetic tape and 
made available for further processing, such as the generation of contour maps. 
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